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We investigate chiral symmetry breaking in the (2+l)-dimensional Thirring model as a function 
of the coupling as well as the Dirac flavor number Nf with the aid of the functional renormalization 
group. For small enough flavor number Nf < Nf 1 , the model exhibits a chiral quantum phase 
transition for sufficiently large coupling. We compute the critical exponents of this second order 
transition as well as the fermionic and bosonic mass spectrum inside the broken phase within a 
next-to-leading order derivative expansion. We also determine the quantum critical behavior of the 
many-flavor transition which arises due to a competition between vector and chiral-scalar channel 
and which is of second order as well. Due to the problem of competing channels, our results rely 
crucially on the RG technique of dynamical bosonization. For the critical flavor number, we find 
iV f cr ~ 5.1 with an estimated systematic error of approximately one flavor. 



I. INTRODUCTION AND SUMMARY 

Introduction. Since the advent of graphene PQ there 
has been an enormous amount of renewed interest in 
(2+l)-dimensional relativistic fermion systems. Various 
variants of (2+l)d quantum electrodynamics [2"HT5] and 
the (2+l)d Thirring model [TM3"3"] are actively discussed 
as effective descriptions of graphene's exceptional elec- 
tronic properties [MKT] . Moreover, they are likewise 
intrinsically interesting: in these models the number of 
fermion flavors Nf serves as a control parameter for a 
quantum phase transition at a critical value 7V f cr . Sev- 
eral previous works provide a substantial amount of evi- 
dence that chiral symmetry breaking (xSB) may be pro- 
hibited even for arbitrarily large coupling if Nf > Nf 1 ' 
[22), I24143T) . This is a similarity to many- flavor non- 
abelian gauge theories in four dimensions which are used 
for particle physics models for dynamical electroweak 
symmetry breaking [421144] . However, the search for the 
quantum critical point has so far been rather challenging: 
in the Thirring model different approximate solutions to 
the Dyson-Schwingcr equations (DSE) yielded values be- 
tween Nf r ~ 3.24 [22] and Nf 1 = oo [23]. By constructing 
an effective potential for the chiral order parameter, up 
to leading order of the 1/Nf expansion iV f cr = 2 has been 
found [26] , Extensive Monte Carlo simulations point to 
iV f cr ~ 6.6 [31 . The true value of _/V f cr is of decisive rel- 
evance for the applicability of such relativistic fermion 
effective theories for condensed matter systems, where 
the number of flavors is typically Nf = 2. If Nf = 2 is 
near iV f cr , the properties of the quantum critical regime 
could take an important influence on physical effects cor- 
responding to dynamical mass generation in the effec- 
tive (strongly-coupled) theories. Understanding such a 
semimetal-Mott insulator transition is also desirable from 
a technological point of view, e.g., with the ultimate aim 
of engineering the band gap in graphene-like systems. 
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In addition to the significant quantitative discrepan- 
cies observed in the literature, a more detailed compari- 
son of the critical behavior close to iV f cr reported in those 
studies reveals our insufficient understanding of fermionic 
field theories in the nonperturbative domain: Kondo [26] 
reports a second-order phase transition with the usual 
power-law behavior as a function of the control param- 
eter Nf. By contrast, in the DSE studies [2U US] an 
essential scaling behavior of the Kosterlitz-Thouless type 
has been found, that is to say, a phase transition of in- 
finite order. As argued in [31] , the nature of the transi- 
tion in these studies appears to depend on whether the 
strong-coupling limit is taken before or after Nf /• 7V f cr . 
The scaling analysis on the lattice [2TJ 132] is consistent 
with a power-law behavior corresponding to a second- 
order phase transition which qualitatively confirms, but 
quantitatively deviates from the Kondo scenario [26] , 

In the present work, we confront the functional renor- 
malization group (RG) approach with the puzzles given 
by the (2+l)d chiral fermion systems at criticality. The 
question is: What is the nature of the xSB phase tran- 
sition being controlled by either the bare coupling g (for 
fixed Nf < N") or the number of flavors Nf (for fixed 
g > g CI )7 What are, if applicable, the corresponding 
exponents determining the critical behavior? Naively, a 
universal answer may not be apparent as (2+l)d fermion 
models are perturbatively nonrenormalizable. However, 
rather than inherent inconsistencies, this difficulty may 
reflect the failure of the perturbative approach [ISMS] , 
In fact, there is a large body of evidence for nonperturba- 
tive renormalizability to all orders in a large- Nf expan- 
sion [T9Tf2"3"] and, more recently, also beyond the 1/Nf 
expansion from a functional RG approach [1H1 ES] ■ As 
a bottom line, (2+l)-dimensional fermion models appear 
to be a paradigm example for asymptotically safe theo- 
ries [5Trf55] that are UV complete due to the existence of 
interacting UV fixed points. 

In the remainder of this section, we briefly summarize 
our findings, starting with a short review of a preceding 
fermionic RG analysis [49]. All technical details and a 
proper embedding and comparison to the literature is 
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FIG. 1. Fermionic RG flow in the Gross-Neveu (GN) - Thirring (Th) coupling plane, compiled from Refs. |49l I50| . The angle 
bisector in the 3rd quadrant defines the theories with pure Nambu-Jona-Lasinio (NJL) interactions. Left panel: for N{ = 1.75 
the Thirring fixed point C lies exactly on the NJL axis with C\ u t oc (S) 2 . Right panel: in the limit Nf — > oo, the Thirring fixed 



point C lies on the Thirring axis with & n t oc (V) 2 . Middle panel: for intermediate flavor number (here N{ 
competition between the scalar-type NJL channel (S) 2 and the vector-type Thirring channel (V) 2 . 



5), we expect a 



deferred to the following sections. 

Competing channels in the (2+l)d Thirring 
model. Conventionally, the Thirring model in 2+1 di- 
mensions is defined by the Lagrangian 

c = c kin + £ lnt - + 2 , (i) 

with Nf flavors of massless four-component Dirac spinors 
ip a , i.e., a — 1, . . . , Nf . The microscopic theory has a chi- 
ral symmetry U(2iVf), cf. Sec.[llJ Depending on the value 
of the four-fermi coupling g and the number of flavors 
Nf, a fermion mass can be dynamically generated, which 
breaks the chiral symmetry spontaneously. Upon inte- 
grating out fluctuations, further interaction terms, being 
compatible with the present U(2iVf ) symmetry, are gener- 
ated by the RG transformations. On the four-fermi level 
there are three further interactions besides the Thirring 
term ~ ("07 M ?/>) 2 : a flavor-singlet pseudo-scalar chan- 
nel (Gross-Neveu interaction), a flavor-nonsinglet scalar 
channel (Nambu-Jona-Lasinio interaction), and a flavor- 
nonsinglet axial channel. However, not all of these terms 
are independent in the pointlike limit: due to the Fierz 
identities we can always choose an arbitrary subset of 
two terms and write the respective other two as a linear 
combination of these. A full basis is therefore given, e.g., 
by the Thirring and the Gross-Neveu interaction. The 
RG flow of the dimensionless couplings in this basis is de- 
picted in Fig.[l] The vertical (horizontal) axis defines the 
theories with pure Thirring (Gross-Neveu) interaction. 
We refer to this axis as "Thirring axis" ("Gross-Neveu 
axis"). The angle bisector in the third quadrant defines 
the theories with pure Nambu-Jona-Lasinio (NJL) inter- 
action, which we refer to as "NJL axis" . In a functional 
RG analysis (the technique is sketched in Sec. |III[ ), we 
find two interacting fixed points A, B, and C besides the 
Gaufiian fixed point O. For the theories defined by the 
microscopic Lagrangian ([I]) being purely on the Thirring 
axis, we find two different phases, which are separated 



by the separatrix through the fixed points B and C (red 
curve in Fig. [I]): If we start the RG flow in the UV with 
a microscopic coupling on the Thirring axis above the 
separatrix, the couplings eventually flow to the nonin- 
teracting Gaufiian fixed point O; if we start below this 
curve, the four-fermi couplings grow large in the IR. In 
the vicinity of the critical coupling, the behavior of the 
system is governed by the interacting fixed point C, as all 
trajectories are initially attracted toward this fixed point. 
We refer to this fixed point as "Thirring fixed point" . We 
associate all trajectories emanating from C with UV com- 
plete fully renormalizcd versions of the (2+l)d Thirring 
model. However, we emphasize that the Thirring fixed 
point only in the strict large- Nf limit lies directly on the 
Thirring axis: even if absent on the microscopic scale, a 
second coupling besides the Thirring coupling will always 
be generated by the fluctuations for finite Nf. 

We now give a heuristic argument for the occurrence 
of the critical flavor number [JHl [SO]- Let us first con- 
sider the (unphysical) flavor number Nf = 1.75, where 
the Thirring fixed point C lies exactly on the NJL axis 
(left panel of Fig. [T]). In this specific case, the NJL axis 
is an IR attractive hyperplane and the RG flow of the 
Thirring model for sufficiently large (negative) coupling 
is dominated by a divergent NJL coupling. A divergent 
scalar-type four-fermi coupling at finite RG scale signals 
bound-state formation in that channel and can be associ- 
ated with the dynamical generation of the corresponding 
order parameter, in this case the chiral order parameter 
(ipij))- We therefore predict that spontaneous breaking 
of chiral symmetry occurs due to a dominance of the 
scalar-type NJL channel. Since the fixed-point positions 
change only smoothly with Nf, we expect this conclu- 
sion to hold also for small deviations from Nf = 1.75. 
By contrast, for very large flavor number, the IR at- 
tractive hyperplane with the Thirring fixed point C is 
located on the Thirring axis, i.e., at vanishing NJL cou- 
pling (right panel of Fig.[T]). In other words, the RG flow 
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FIG. 2. Log-plot of the dynamically generated renormalized 
fermion mass mj^ f in units of the renormalized mass mj^ „ of 
a radial bosonic excitation on top of the chiral condensate. 
Our estimates for the critical flavor number, N" ~ 5.1 for a 
linear regulator is amended by an error bar (gray shaded area) 
obtained by varying the regularization scheme (e.g., N" ~ 5.8 
for the sharp cutoff). 



is governed by a strong vector-channel dominance, which 
generically inhibits x^B. For intermediate flavor number 
(middle panel of Fig. [lj , we thus expect a transition be- 
tween the scalar- type NJL channel, which triggers xSB, 
and the vector- type Thirring channel, inhibiting x^B. 

Dynamical boson fields. A more quantitative pic- 
ture of this quantum phase transition is difficult to obtain 
in the purely fermionic language with point-like interac- 
tions. A quantitative RG analysis requires the inclusion 
of dynamical chiral and vector bosonic degrees of free- 
dom, see Sees. |IV|VI[ in order to study the interplay of 
these competing channels as a function of Nf. This is the 
objective of this paper. In particular, we will show how 
the competition between the NJL (dominant for small 
Nf) and the Thirring channel (dominant for large Nf) 
leads to a decrease of the dynamically generated fermion 
mass with Nf, culminating in a complete vanishing at 
iV f cr ~ 5.1 (linear regulator); cf. Fig. [2} We can also map 
out the order of the phase transition as a function of 
Nf'. in Fig. [3] we depict the order parameter {tp) oc (f/>V) 
vs. N" /Nf — 1 in a double-log plot, showing very good 
compatibility with a second order phase transition with 
scaling behavior {<p) oc (Nf r / Nf — l) b and the univer- 
sal critical exponent b ~ 0.44. To avoid any possibility 
of confusion, let us denote the critical exponents for the 
phase transition as a function of Nf (at fixed overcrit- 
ical coupling) with Latin letters, and the ones for the 
phase transition as a function of the coupling (at fixed 
Nf < iV f cr ) with Greek letters. We find that the lat- 
ter phase transition is also of second order, the critical 
behavior of which is determined by the RG flow in the 
vicinity of the Thirring fixed point. We give our predic- 
tions for the universal critical exponents v, uj, and for 
various flavor numbers in Table [IT] on page [15] 

In addition to our quantitative results, it is one of 
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FIG. 3. Order parameter (black dots) vs. N"/N f — 1 in a 
double-log plot, showing very good compatibility with power- 
law scaling behavior {tp) oc (Nf T /N{ — l) b . The slope of the 
regression line (gray) is b ~ 0.44. 

our most important observations that a simple partial 
bosonization a la Hubbard-Stratonovich is not sufficient 
for resolving the competition between the various chan- 
nels. A method to deal with this challenging problem is 
dynamical bosonization, as detailed in Sec. |VII| As there 
is already a large variety of analytical as well as numeri- 
cal studies on the (2+l)d Thirring model in the literature 
- with partly contradictory results - we finally perform 
a careful and critical comparison between the literature 
and our work in Sec. |Xj 

II. CONDENSATION CHANNELS AND FIERZ 
BASIS 

Let us start by presenting our conventions for general 
fermionic models in d = 2 + 1 Euclidean dimensions with 
local quartic self-interaction, being invariant under the 
full symmetry of the noninteracting theory. We con- 
centrate on representations of the Dirac algebra as they 
occur naturally in the language of effective theories de- 
scribing electronic interactions on graphene's honeycomb 
lattice [3~4"H4"0] and in cuprates 56 -f6"T] . 

The models shall satisfy Osterwalder-Schrader positiv- 
ity [62], requiring invariance of the action under gener- 
alized complex conjugation defined by tp' := i-07 3 and a 
simultaneous reflection of the Euclidean time coordinate, 
which we choose to be X3. For a detailed discussion of 
our chiral conventions, see Refs. 63, 64J. Consider the 
4x4 reducible representation 

7p=(^T"), = 1,2,3, (2) 

of the Dirac algebra {7^,71/} = 2<5 MZ ,, reminiscent of 
the chiral representation in four dimensions. Here, 
{c A i} M =i,2,3 denote the standard 2x2 Pauli matrices. In 
this formulation, the Dirac fermions thus have four com- 
ponents. There are now two other 4x4 matrices, which 
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anticommute with all 7^ as well as with each other, 
1 2 



' h 

Together with 



and 75 = 71727374 



la 

-1 2 



(3) 



I4, 7^ : = 2^'^"] 0" < i > 74 ' i7 ^ 75 ' 745 := i7475 ' 

(4) 

these 16 matrices form a complete basis of the 4x4 Dirac 
algebra, 

{ia} a=1 i 6 = {I4, 7^,74, 7^- 74, 17^75, 745, 75} ■ 

(5) 

The general fermionic Lagrangian compatible with 
U(2iVf) chiral as well as C, V, and T discrete symme- 
try has the form QUI HQ [Ml \M EH] 



C = V a i« 



91 



93 

2N< 



.94 

2N f 



{Af 



(6) 



with the bare couplings 31,.. ..4 carrying an inverse mass 
dimension, and with the flavor singlet channels 

(vf : = (r^r) 2 , (p) 2 ■■= (^w a ) 2 , (7) 

and the flavor nonsinglet channels 

(S) 2 := (^ b ) 2 - (V>74^ b ) 2 - (^W) 2 

+ (^<W) 2 , (8) 



(9) 



Here, we have abbreviated (ip a ip b ) 2 = ifj a ifj b ijj b il! a , etc. 
JVjf denotes the number of four-component (reducible) 
Dirac spinors, such that a, b = 1, . . . , N{. By means of 
the Ficrz identities, 

{Vf + (S) 2 + (P) 2 = 0, -4(V) 2 - 3(5) 2 + (A) 2 = 0, 

(10) 

any two four-fermi terms can be rewritten as a linear 
combination of the remaining two. Put differently, by 



adding a linear combination of Eqs. ( 10 ) with coefficients 



a, G R (in units of some inverse mass scale) to the La- 
grangian (|6|, 

£ = + JLfo + ai - Aa 2 )(V) 2 

+ ^(.92 + a x - 3a 2 )(S) 2 + ^(§3 + «i)(^) 2 



+ ^(54 + « 2 )(^) 2 , 



(11) 



the a,i are redundant parameters: in a full computation 
of the functional integral, any physical quantity has to 
be independent of aq. This no longer necessarily remains 
true once approximations are employed. A particular 
example is given by mean-field theory, where this so- 
called "Fierz ambiguity" or "mean-field ambiguity" can 
have a sizable influence on the results, limiting its quan- 
titative reliability |65| . The ambiguity is absent in the 
purely fermionic renormalization group equations with 
momentum-independent couplings |49j . A solution of the 
Fierz ambiguity using the functional RG in a partially 
bosonized setting [66] requires dynamical bosonization 
[6 7) as will become important below. An alternative ap- 
proach in the purely fermionic description employs a new 
parametrization of the momentum structure of the four- 
fermi couplings, see 

The particular choice of couplings g and g as used in 
[49] is recovered for ax = —g 2 — 3<?4 and a 2 = —34 and 
the definition 

9 : = .9i ~92 +34, (12) 
9 ■= -92 +.93 ~ 3.94- (13) 

The Thirring (vertical) axis in Fig. n] corresponds to g, 
whereas the Gross-Neveu (horizontal) axis is associated 
with g (more precisely with their dimensionless counter- 
parts). Upon choosing ax = —33 and a 2 — —.94, the 
Lagrangian reads 



- (r^ b ) 2 - (?W) 2 + (^w 

where we have defined the new couplings 



9v ■= 

9<$> - 



-.9i - 
2(32 



33 - 
- 33 



4 34 = 9 - 3, 
- 334) = -2g. 



(14) 

(15) 
(16) 



This form is convenient in order to investigate the compe- 
tition between the vector (V) 2 and NJL-type (S) 2 chan- 
nel for Nf > 2. For Nf — 1 one might however choose yet 
another basis, 



C 



Nf 



2.92 - .93 



231 - 93 



334 



334 



(ipip) 2 - (^74-0) 2 - {ipi5ipy 



(17) 

The Dirac spinors can be projected onto their Weyl com- 
ponents, using the projectors^] 



P 



(45) 



L/R 



2(1 ±745) 



l a (-i)la 
il 2 1 2 



(18) 



1 We note that one could choose more adapted representations in 
which 745 = diag(l2, — 12) and thus the Weyl spinor \ a (x a+JVf ) 
is simply given by the upper (lower) two components of the Dirac 
spinor ij) a . The "graphene representation" |35j|36| is of this type. 
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where the last equation holds for the representation 
([3]). This suggests the decomposition [TO] 



,a+N { 



v ~ V2 V(x a -x a+Ni ] 



(19) 



Tjj a = — (x a -x a+Ni , H) (X a + X a+Nt )) , (20) 

a = 1, . . . , N[ , with the definition of the Dirac adjoint x 
chosen such that x^ := l X (J 'i m agreement with ip^ := 
iip"t3- I n other words, we can trade the N{ flavors of four- 
component spinors ip for 2Nf flavors of two-component 
spinors \. Therewith, the theory (14) can equivalently 



be described using an irreducible representation by an 
action consisting of 2Nf flavors of two-component Weyl 
spinors \, X, 



^(xV^) 2 + H(xV) 2 , (21) 



where we have introduced the collective indices i, j, run- 
ning over 2N{ flavors, i,j = 1, . . . , 2N{. It is this repre- 
sentation in which the U(2JVf) symmetry is manifest, 

U(2A f ): x l ^U t] X \ X z ^X°(uy\ U G U(2iVf). 

(22) 

Below, we use this formulation, allowing us to conve- 
niently introduce collective low-energy degrees of free- 
dom. 



III. FUNCTIONAL RG APPROACH 

The functional RG has become a standard method to 
investigate strongly-interacting field theories. In partic- 
ular for critical phenomena, it is a versatile tool yielding 
quantitatively accurate results in many cases. Conceptu- 
ally, the functional RG can be formulated in terms of an 
RG flow equation for a generating functional. Among 
the various formulations, the Wetterich equation [69] 
representing the flow equation for the effective average 
action has become the most widely used method ow- 
ing to its flexibility, numerical stability and direct ap- 
plicability to physics problems. For reviews see [70Tf76] 
and [TTHT9"] for a particular emphasis on fermionic sys- 
tems. Prominent benchmark examples in three dimen- 
sions are bosonic O(N) models [701 [HOI EI] or the Gross- 
Neveu model [55l E21 E3] . 

The effective average action is a scale-dependent 
variant of the standard generating functional for 1PI cor- 
relation functions. The scale k denotes an IR regulator 
scale separating the UV modes with momenta p 2 > k 2 
which have already been integrated out from the IR 
modes with momenta p 2 < k 2 which still have to be aver- 
aged over. Tfc is constructed in such a way that it can be 
related to the microscopic bare action — > 5a if k 
approaches the UV cutoff scale A, while it approaches the 
standard full effective action in the IR limit —> T. 



The effective average action obeys the Wetterich equa- 
tion [69] . 



d t T k [<f>} = -STr 



d t Rk 



t = *A, ( 23 ) 



where the trace runs over all internal degrees of free- 
dom (flavor, spinor, momentum) as well as field degrees 
of freedom. Here, the field $ represents a collective 
field variable including all bosonic or fermionic fields un- 
der consideration. The denominator contains the second 
(2) 

functional derivative T k [$] with respect to the field <& 
together with the regulator function Rk which can be 
thought of as a momentum-dependent mass term. 

Once the initial condition of the flow is fixed in terms of 



the microscopic action 5a, the exact solution of Eq. (23) 



provides us with an RG trajectory of in the theory 
space of action functionals, the end-point at k = of 
which is the exact effective action T. 

As it is, in practice, difficult to find exact solutions, 
the Wetterich equation can also be used to find approx- 
imate solutions by means of systematic and consistent 
expansion schemes of the effective action. While per- 
turbation theory constitutes one such expansion scheme, 
nonperturbative schemes based on vertex or operator ex- 
pansions are equally legitimate and clearly superior at 
intermediate or even strong coupling. For critical phe- 
nomena and the analysis of long-range order, derivative 
expansions in terms of the order-parameter fields have be- 
come a standard tool yielding accurate results in many 
nontrivial examples. In the present work, we will also 
use a truncation of the effective action in the spirit of the 
derivative expansion. However, as the chiral order pa- 
rameter field is bosonic, the derivative expansion has to 
be set up not only for the microscopic fermionic fields as 
in [49 , but requires collective bosonic field variables, as 
introduced in the next section. As a side remark we note 
that continuation of the RG flow into regimes with broken 
symmetry can also be achieved in purely fermionic de- 
scriptions without Hubbard-Stratonovich transformation 
by inserting an infinitesimally small symmetry-breaking 
component in the initial action [84j [85] . 

The above mentioned regulator function Rk is a to 
some extent arbitrary function, satisfying a few con- 
ditions for implementing a meaningful regularization 
[701 ESI E2- Whereas exact solutions do not at all depend 
on the specific choice of Rk, approximate solutions and 
even estimates of universal (i.e., regularization scheme in- 
dependent quantities) can depend on Rk- In the present 
work, we will use this variation of universal quantities as 
a function of Rk as an estimate for the systematic error 
introduced by our approximations. We consider results 
obtained for the linear regulator as our best estimate, as 
it satisfies RG optimization criteria [73, 86]. By contrast, 
the sharp cutoff is known to introduce strong regulator 
artifacts which we therefore use to maximally span the 
error bar for our approximations. 
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IV. LOW-ENERGY DEGREES OF FREEDOM 



1, . . . , iVf, with the breaking pattern 



The partition function of the theory defined by (21 ) 
Z = J VxV X exp(-S), 



(24) 



is in fact equivalent to the partition function of the 
"mesonic" theory of a (2iVf) x (2iVf) scalar matrix field 
(f)^ and a vector field V^, coupling via a Yukawa- type 
interaction to the fermions, 



Z = N J VxDxDcfiDV exp 



(25) 



where i,j 



1, 



, 2N[. The equivalence can be 



seen by multiplying Eq. ( 24 1 with appropriate GauBian 
(Hubbard-Stratonovich) factors, 

1 = TV J V<j> exp ' 



( _ ,44 .h 



1 / _ „ hv - t , 
■- myV^ - —x o-^x 

2 \ my 



1 = N J VV exp 



(26) 
(27) 



with some normalization constants A/", not affecting any 
expectation values. The scalar matrix field is Hcrmitian, 
4>' = <j>, and the vector field V is real. The four-fermi 
terms in Eq. ( 24 ) are then precisely canceled if the con- 
straints 



K_ 

2m 2 



9±_ 
2N t 



h 2 y 
2m v 



9v_ 
2N{ ' 



(28) 



are imposed at the microscopic scale. From Eqs. (26) 



(27) we can read off the properties of the boson fields 



under chiral transformations, 



5 lJ 

j« Tr. 

2N t 



ip := Tr ( 



(30) 



The ip field is parity odd and can be attributed to the 
parity breaking channel (P) ~ V> a 745?/> a = x^X 1 ■ By con- 
trast, a vacuum expectation value of the traceless part $ 
corresponds to the dynamical breakdown of chiral sym- 
metry, 

± o & (ff) = (x a x a - x a+Nt x a+Nt ) * o, 

(31) 



U(27V f ) ->■ U (JV f ) ® U (A/f) 



(32) 



We can trade the (2Nf) x (2iVf) Hermitian traceless ma- 
trix $ for its independent components $ Q , 



*« = V^^)* 



(33) 



i, j = 1, . . ., 2iVf, a = 1, . . . , (2iVf) 2 — 1, where the t a are 
the generators of SU(2iVf) in the fundamental represen- 
tation, normalized so that Tr(t a tp) = 8 a p/2. 

In the one-flavor case Nf = 1 this formulation is 
equivalent to a partial bosonization of the Fierz basis 
in the four-spinor representation (17): the Hubbard- 
Stratonovich transformation leads to the equivalent 
Yukawa-type theory with three scalar modes (cr, r, tt) ~ 
(iptp, 4>j4i/j, V'75'0) and a vector mode V M ~ tpj/iifi with 
Lagrangian density 



1 



U(2iV f ): 0« ^U tk cf> kl V^V^, (29) J 

[7 G U(2iVf). The scalar matrix </> may be decomposed 
into a traceless part and its trace [ST] 



C Ni=l = ipitfip + ^m 2 (a 2 + t 2 + tt 2 ) 

+ ift, CT V> (f + i74T + i7 5 7r) ^ - JivV^^ip. (34) 

From the discussion of the fermionic RG flow [?9] we ex- 
pect that the long-range dynamics of this system is dom- 
inated by the scalar NJL-type channel. For large enough 
coupling /i 2 /m 2 we thus expect the scalar mode to ac- 
quire a nonvanishing vacuum expectation value (VEV), 
e.g., in the a direction, and the spectrum in the broken 
phase consists of two massless Goldstone modes, e.g., r 
and 7r, and a massive radial mode a. The corresponding 
critical behavior is an interesting problem by itself: this 
system can be viewed as an effective low-energy theory of 
spinless electrons on the honeycomb lattice and thus as 
a simple model for suspended graphene [33]. It has also 
been studied in the context of magnetic catalysis [55J . In 
the following, we will focus on the case Nf > 1 where a 
true competition between the two channels (V) and (5) 
is expected. 

In contrast to the purely fermionic formulation [49) . 
the bosonized formulation presented here is well suitable 
to quantitatively describe the spontaneous breaking of 
chiral symmetry. Loosely speaking, the bosonic fields 
M and V/j, parametrize the possible formation of bound 
states of the fermionic fields ~ x l X^ an d ~ X l o>X% re- 
spectively. The corresponding critical phenomena of such 
a strongly-correlated system require nonperturbative ap- 
proximation schemes. The functional renormalization 
group formulated in terms of the Wetterich equation is 
such an appropriate tool and has already shown its quan- 
titative reliability in other (2+l)-dimensional relativistic 
fermion systems, see e.g., [5511521155] . In the effective ac- 
tion we then have to take into account also higher boson- 
boson interactions generated through fluctuations, e.g., 



(Tr$ 2 ) - ($ Q 1> Q ) 2 , 
1 



(35) 



Tr$ 4 = 



2Nf 



+ 2d af 3 e d 7 6e<f> a $p% ( Z , 8, (36) 
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where the d a p^s are the structure constants for the group 
SU(2A/f). For Nf = 1 the d a ^ 7 's vanish. For computa- 



tional reasons, we use the representation ( 25 ) for the case 
of general Nf with scalar matrix field (j? 3 having arbitrary 
trace. 

We parameterize the dynamics of the low-energy de- 
grees of freedom by an ansatz for the effective average 
action T^. This ansatz corresponds to a systematic ex- 
pansion of the action in powers of the field gradients to 
second order. Moreover, we simplify the discussion by 
considering interactions only up to fourth order in the 
fields. Our ansatz for the effective action then reads 



T k = 



J V,k 



2 

A v .k 



2 -VM 



®— (yV y 1 ) 2 + V— (yV'V 



(37) 



where := <9 M K - dyV^ and i,j = l,...,2N{. 

Uk(4>) describes an effective potential in the scalar sec- 
tor. All couplings in the effective action are understood 
to be scale-dependent, indicated by the index fc. It is 
straightforward to match this action exactly to the micro- 
scopic models discussed above. For instance, by setting 
all bosonic couplings to zero at the UV cutoff scale k = A, 
we return to the fermionic action (21). Alternatively, we 
can set the fermionic couplings to zero, gv,A,9(j>,A = 0, 
and choose the compositcness conditions 



Z Vtk ,A V u, ( k -> 0, for k -t A, 



(38) 



which guarantee that the bosonic fields have no kinetic 
term and are purely auxiliary at the high scale. Fur- 
thermore, satisfying the constraints ( 28 ) and setting the 



higher-order bosonic interactions to zero at k — > A leads 
us to the Yukawa-type theory as given in Eq. (25 1. 



This implements the Hubbard-Stratonovich transforma- 
tion at a fixed scale k = A in the functional RG con- 
text. Of course, as soon as we start integrating out 
modes, the couplings set to zero at one scale can be 
generated at lower scales due to fluctuations. Also the 
constraints (28) will generally not be satisfied at lower 
scales. Whereas the purely fermionic flow has been ex- 
tensively discussed in |49j . we will now concentrate on 
the partially bosonized variants. We start with comput- 
ing the RG flow in Sec. |VI| for the fixed- field variables in- 
troduced above. As the Hubbard-Stratonivich transfor- 
mation demonstrates, the effective action (37) contains 



some redundancy, as field reparametrizations can trade 
various operators for one another. In Scc. |VIll we exploit 
this redundancy to introduce scale-dependent fields that 
adjust the fermion-boson couplings dynamically in order 



to effectively perform Hubbard-Stratonovich transforma- 
tions on each scale (dynamical bosonization) . In the fol- 
lowing, we will mostly ignore the momentum-dependent 
terms in the vector channel, i.e., we will consider the 
pointlikc limit 

Z v .k -> 0, A Vtk -> 0, Cfe 0, for all Jfc. (39) 

The beta functions will however be computed for general 
Z v ,k,A v ,h > (but Ck = 0). 



V. SCALAR MASS SPECTRUM 

Due to U(27Vf) symmetry, the effective potential Uk(4>) 
necessarily has to be a pure function of U(2iVf )-invariant 
quantities. The Hermitian scalar matrix field can be di- 
agonalized by a U(2Af) rotation [cf. Eq. (29)], 



4> i — ^ [701^ 



(m x 
m 2 

\0 



\ 




m 2 N t J 



(40) 



with real eigenvalues fhi. In a quartic approximation 
Uk(4>) can then be parametrized in terms of the two in- 
variants p and r, 



p = — Tr <f) 2 = 



t = — Tr 



2N f 



E 



^4 

m 4 — 



AN f 



(41) 

2 

(42) 



At higher field orders the potential can depend on addi- 
tional (suitably defined [89 ) higher order invariants con- 
structed from f„ - Tr(0 2 /2 - p/2N f ) n for n > 3. We 
expand the scalar potential about its fc-dependent min- 
imum (po,fc7 T o.fe)- In the symmetric (SYM) regime the 
potential is minimal at the origin (po,k, T o,k) = (0,0), 
whereas in the chiral symmetry broken (%SB) regime we 
allow for a nonvanishing VEV which we assume to be ac- 
quired along the p direction, i.e., po,fc > an d tq k — 0: 



Uk{p,r) 



m l,kP- 



^rP 2 + A 2 ,feT, SYM regime, 



HjMp - Po,k f + A 2: feT, xSB regime. 



(43) 



If the potential solely depends on p and r, it is sufficient 
to evaluate its flow equation in a two-dimensional sub- 
space of all possible scalar configurations. Wc consider 
the traceless class, valid for Nf > 1, 



b 13 ) =: mdiag(e, +l ,+1 ,-e, 

Nf — 1 times 



-i, 



■1), (44) 



Nf—1 times 
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TABLE I. Spectrum of the scalar mass matrix 5 2 U k /6(j)54> for d 2 U k /dr 2 = d 2 U k /dpdr — and Nf > 1. Primes denote partial 



derivatives with respect to p, L/^ = d n U k /dp n : and U k}T = dU k /dr. For fh = fh(p, t) and e = e(p, r) see Eqs. (|45|)-(|46 1 



Eigenvalue 



Degeneracy 



>u t -2-Uk, T |3e 



1 + 



u 'k + i^|47V f (7V f - i)U£ + (4 - N f )U k , T + [4N { U k ' + (3iV, - 4)t/ fc , T ] t 

±[{4N { (N l - 1)E#' + (N t + 2)U k , T ) 2 

+2(l6Y f 2 (Y f - l)U' k ' 2 - 4N f (N t + 2)(3JV f - 2)[/£t4, T 

-(3Y f 2 - 42V f + 4)[/ fe 2 T )e 2 + (4AT f [/£ + (3N ( - 2)U k , T f e 4 ] ^ J 



2 



fc T — Uk,r 



2 Ufc . T JV. 



1 

1 + 1 



4(7V f - 1) + 4(iV f - 1) 
27V f 2 — 4Nf + 1 
2 (7V f - l) 2 



where e > 1, m € R. For any p and r with r/p 2 
(Nf — l)/4JVf, the parameters m and e are given by 



< 




e 2 = l 



1 



4JV f 



- 1 



(45) 



(46) 



iVf — 1 p 2 



If the flow eventually chooses a vacuum configuration 
with e = 1 and fh > 0, the chiral symmetry is sponta- 
neously broken while parity symmetry remains preserved, 



Po ft 



i') = f IS -VI'"' 



» G*X* - x' +lv 'x° + '*'> = 7 s 0. (47) 



For simplicity, we assume in the following that 
d 2 Uk/dr 2 = d 2 Uk/dpdr — 0, which holds in the case 



of the quartic approximation ( 43 1 . The spectrum of the 
scalar mass matrix S 2 Uk/8cj>5(j) for Nf > 1 is given in 
Table D 

For a vacuum configuration with e = 1 we have To,jfc = 
0. In the SYM regime po,fc = 0, all modes are degener- 
ate and have mass = 9t/fe/9/3|( /J , T )=(o,o)- I n the xSB 
regime with p a . k > and dU k /dp\ {p ^ = ( po fc:0) = we 
find 2N{ 2 massless scalar modes corresponding exactly to 
the number of broken generators in the symmetry break- 
ing pattern 



be a general result, holding also beyond our quartic ap- 
proximation ( 43 1 , as long as no higher-order invariants 



t„ become important. 



VI. PARTIALLY BOSONIZED RG FLOW 

Next, we determine the flow of the model in a par- 
tially bosonized description. This corresponds to start 
at the high scale with the Hubbard-Stratonovich trans- 
formed action of ( 25 1 and the compositeness condition 



( 38 1 . As is standard in this type of truncation, higher- 
order fermionic interactions are set to zero also on all 
lower scales. 

We fix the standard RG invariance of field rescalings 
by defining the renormalized fields as 



V„ 



- 7 1 / 2 



W,k "> 



Vu, 



<A — Z 1 / 2 -^ 

X — \ t x 



(49) 
(50) 



U(2Y f 


-> U (Nf) <g> U (N t ) , 


(48) 


Av.k '— Z v 1 Ay t k, 


V 


'— Z 4> 1 k Z VM kd 


-4- 
Vk, 


in accordance with 


Goldstone's theorem. 


Addition- 


2 . ry-l 1—2 - 2 

m v — ^v,k K m V,ki 


hi 


'■= Z 4>\ Z x% kd ~ 


4 h 2 


ally, we obtain one 


massive radial mode 


with fh 2 = 


p := Zy 2 k k d - 4 p k , 


h 2 v 




n V,k- 



The dimensionless effective potential for space-time di- 
mension d then reads 

u(p,f) := k- d U k (Z^ k k d - 2 p,Z^ 2 k k 2 ^r), (51) 

where p := Z<^^k 2 ~ d p and f := Z 2 k k 2{ - 2 ^ d ^T, and the 
dimensionless renormalized couplings are 



(52) 
(53) 
(54) 



2po.fc<9 Uk/dp and 2Nf — 1 massive modes in r direc- 
tion with to 2 = (p 0tk /Nf)dUk/dT. Since St/ =0 By evaluating the Wetterich equation p3) for a constant 
we expect the 2iVf 2 — 1 degeneracy of the masses to 2 to scalar background field (44 1, we obtain the flow of the 
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dimensionless scalar potential 

d t u = —du + (d — 2 + r)(j,)pu' + (2d — 4 + 2j]^)fu t f 

+ 2v d ^V,; m) + 2v d d£ { * )d (m 2 v + vp- r, v ) 

»=1 

-4^d 7 [(iV f -l)^ )d (m 2 ^ ;7?x ) 
+ l^ d (m 2 e 2 h 2 ;r lx ) 



with u — u(p,f), v! = du/dp, 



m 



2._ 



2 " d TO 2 



(55) 

= du/df and 
The dimensionless scalar masses 



mf = m 2 (p,f) can straightforwardly be deduced from 
Table HI We have further defined the anomalous dimen- 
sions 



V<p/ x /v 



-d t In Za 



'x/V,k- 



(56) 



The flow involves the threshold functions £q B ^ F ' (...), 
which encode the details of the regularization scheme; 
their definitions and explicit forms for the linear and the 
sharp cutoff are listed in the Appendix. We have ab- 
breviated v' 1 = 2 d+1 n d / 2 T{d/2), i.e., Ug 1 = 8tt 2 . As 
discussed above, we work here with the two-component 
Weyl spinors such that the dimension of the gamma ma- 
trices is d 7 = 2. By suitable differentiation of Eq. (55) 



we obtain the flow of the scalar couplings occurring in 



Eq. (43). In the SYM regime, 



d t m\ = d t v! 
d t \i = d t u" 



(p,f) = (0,0) 



(p,f)=(0,0) 

whereas in the x^B regime, 
1 



<9 f A 2 = d t ufi 



(p,f)=(0,0) 



d t K = -^-d t u' 



(P,t) = (k,0) 



d t Xi = d t u'' 



(p,f) = (K,0) 



d t X 2 = d t u : f 



(57) 
, (58) 

(59) 



(p,f)=(«,0) 

(60) 



with the dimensionless VEV k := Z < j ) j : k 2 ^ d pQj £ and the 
abbreviations vy 1 ' = d"u/dp n and = du/df. 

Similarly, the flow equations for all other couplings 
present in the effective action are straightforwardly ob- 
tained by suitable projections of the Wetterich equation 
( [23] ) . This amounts to a summation of all possible f-loop 
diagrams where the vertices are given by full (though 
truncated) vertex functions and the inner lines corre- 
spond to the full propagators. With the useful Math- 
ematica package DoFun [§U], the evaluation of the RG 
flow equation can be automated easily. The beta func- 
tion for the vector mass reads 



d t my = (-2 + r) V )my - 2v d (d + 2) ^ )d (m v + vk\ r\y)p 



2v d 



2N { 2 4 B)d (u'; m) + (2N { 2 - l) tf )d (u> + f« f ;^) + if )d {v! + 2nu"; ^) 

. 7 ^ (i^h^rj^h^hy, 



8v d {d - 2)d 1 N l (F)d , 2 
-j n \w t n <t» T lx) n v 



(61) 



with the derivatives of the potential u^ n ' and u t ? being evaluated at the minimum (p, f) = (k, 0). In the symmetric 
regime we have of course k — 0. For the vector- vector interaction p and the vector-scalar interaction v we get 



(A 1X1 ^ 1 2^(rf 2 + 10d+12) (B)d 2 2 
o t p = (d — 4 + 2rjy)p H TT~H ^2 \ m v + VK 'i VvjV 



d + 2 



2v d 



2N { 2 4 B)d («'; tit) + {2N { 2 ~ 1) lf )d (u' + fiuMt) + ef )d (u' + 2ku"; V<t> ) 



JB)d 



(B)d, , 



16v d {d-2)(4-d)d 1 N i {¥)d K 2 s ,4 128v d {4-d)d 1 (F )d ( K , 2 U 2.4 

d(d + 2) 2 ^V%)V+ ^(^2) 3 m^w^ 

d(d + 2)N { K h UV%)W, 



(62) 
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d t v = (d - 4 + 77^ + r\v)v + 2(d + 2)v d £^ d (my + vk; r}y)[iv 



+ v d 



(2N 2 + 1 ,, 2N f 2 - 1 

U H TTT^j U.f 



V ^f 2 



27V f 



2N i 2 if )d {u';^) + {2N { 2 - l) tf )d (u' + 



+ 4 B)d ( U ' + 2^";77 ) 



1 



+ 



2vd 
N f 2 



m 2 e (BB)d( u ^ m 2 y + VK . ^ w) + (2 N 2 _ y ^(BB)d (u , + ^ + ^. ^ ^ 



+ ^^^(ti' + 2ku", to^ + ^k; ?70, 



2 1 ^V^U^*) W 



dN f 



lKg>\%hl; V X ) h 2 v hl - 32Vd [ 4 \ T2 d)d ^ 2 4P*&hl; Vx ) h v h% 



dNf' 



(63) 



In the symmetric regime the flow of the Yukawa coupling h<j, is unambiguous. In the broken regime however the 
diverse scalar modes in general can develop different couplings. We will focus on the Goldstone-mode coupling to the 
fcrmions, which is expected to give the dominant contribution for aspects of criticality. However, this Goldstone-mode 
projection may introduce artifacts deeply in the broken regime, see below. The flow equation reads 

d t h 2 = (d - 4 + jfo + 2 Vx )hl 



+ 



N f 

H ttk 

N f 2 



ti 



47 )d (^« , + ^«,f;^.^)-4T )d (^^,« , +w , ;% ,fj ) 

r 2m 2 u"^ )d ( f-h%u'- Vx ^) + ((2Nf 2 l)u f + 2(N f l)u") £^ d (f-h 2 ,u',u> + %u f ;V x ,V4>) 
+ (2u'' + u f )ll^ )d {^,u',u' + 2^''^^)] h%-8v d d£^ )d {^h%m 2 v + VK-,in x ,m)hlh 2 v . 



(64) 



For the flow of the fermion- vector coupling we find 
dth 2 , = {d-A + r] V + 2r] x )h v 



4v d (d - 2) 
dNi 



2N f 2 £^ )d (^h 2 ,u';r lx ,^) + (2N { 2 - l)£^\f-h%u> + ^u f ; Vx , V4> ) 



(FB)d / K ,2 / , K 



+ e { ™ )d (jt { h%u' + 2 K u";r 1x ,r 1<i ,) 



2.2 



hih 



8v d 
dN f 2 



2N l 2 £^ )d (^h 2 ,u / ;r lx ,r l4> ) + (2Nf 2 - l) 6™ )d { f-h%v! + f-u, f ; Vx , V<j> ) 

+ ^f )d (^h 2 ,u' + 2 K u"; Vx ,^) 
8v d (d~2) 2 JFB)d, K ,2 2 



Kh v 



d vjv ' 

and the anomalous dimensions read 



H,i {f-h^my + vn;r, x ,ri v )h v - 



4 _lQv d {d 2)^(FB) d( ^ /i 2 )m 2 / + ^ ;??x)r?y)W) (65) 



dNf 



+ [Nf 2 u" 2 + (Nf - l)u. f 2 ] mfl d {u\ v! + 2ku"; V<j> ), 
dNf 

4vd 
dNf 



Vx = 



2Nf 2 rnf^ d {f-h%u'-,r, x ,^) + (2N { 2 - l) m£f )d (^/4 , u' + f-u , f ; Vx , ^) 



+ m{ i,2 )d (w { h l> u ' + 2ku "> ^) 



K\ - 8w d m 1 F 2 B)d ( ] ^/i|,my + vn\j] x ,j] V ) h' v , 



Vv 



I6v d (d 2)d lN f mf)rf( hl + l6vdd ^ m{ p d{ « h l,r, x )hlh 2 v . 



(66) 

(67) 
(68) 



Within our truncation we find that the flows of the two possible vector-field kinetic terms are in fact equivalent, 
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d t A k = d t Z 



agonal, G 



(V) _ 



v,k- The vector propagator thus is di- 
liVi k — ^^/( z v,kP 2 + ™v.k)- The definitions of 
the threshold functions e {B/F)d (. . ■ ) and m ( : B / F *> d (. . . ) are 
listed in the Appendix, together with their explicit forms 
for linear and sharp cutoff. The flow equations have been 
independently verified for the symmetric regime in local 
potential approximation with the DoFun package |90) . 

As is the case in the fermionic formulation, we expect 
the flow to be dominated by the vector channel for suf- 
ficiently large flavor number Nf. Whether for large cou- 
pling the V field then can develop a finite vacuum expec- 
tation value is an interesting question on its own right: 
e.g., (2+l)-dimensional models exhibiting spontaneous 
breaking of Lorentz symmetry have been investigated in 
[9"1"H9"5] . By contrast, one could suspect that the vector 
mass eventually flows to zero, triggering a close resem- 
blance of the strongly coupled Thirring model for large 
flavor number to a U(l) gauge theory. This is in fact 
the prediction of the large- Nf studies [3TJ [SB]. In the 
present work, we are mainly interested in a possible scalar 
condensation corresponding to chiral symmetry breaking 
and thus leave this issue for future studies. 

As long as the vector mass m v does not become too 
small it is then sufficient to consider the pointlikc approx- 
imation in the vector channel, Zy,k ~> 0, corresponding 
to my — > co. In this limit, the beta functions no longer 
depend on h v and m v separately but only on the ratio 
gv = Nfh v /m v , reflecting the RG invariance of fields 
rescalings. In the remaining flow equations the vector 
anomalous dimension r\y completely drops out. The flow 
of Z Xt k given by Eq. (67) is driven by a competition be- 
tween a positive scalar loop term cx h\ and a negative 
vector loop term oc h v . In the pointlike approximation 
of the vector sector, m v — > co and hence the vector loop 
term vanishes. For reasons of consistency, we will there- 
fore suppress also the scalar loop term in Eq. ( |67| , i.e., 
we treat the fermionic sector in the leading-order deriva- 
tive approximation rj x = 0. This is compatible with the 
observation that the flow of the fermionic wave function 
renormalization in 3d fermion systems at criticality is 
usually very small [SSJ [El [Ml [S7HSS]. As we shall see 
in Sec. |VII[ this assumption is exactly fulfilled for large 
number of fermion flavors Nf — ¥ co, where the fixed-point 
equations can be solved analytically; in fact, this is the 
known result of the 1 /Nf expansion 20J . 

In comparison, the scalar anomalous dimension is non- 
vanishing in this limit, and we thus expect the flow of 
Z^f. to be crucial also for finite Nf. We note that sim- 
ilar observations have also been made in other Yukawa- 
type systems in three dimensions, where a nonvanishing 
?70 ~ 0(1) is essential in order to find the correct critical 
behavior [SFJ [HI [Ml W\ ■ Also the recent results for the 
Nf = 1 model [33] fit into this scheme. 

For the remainder of this section, we omit the vector- 
vector interaction fi and the vector-scalar interaction 
v for simplicity. We concentrate on the UV structure 
only in order to compare it with our previous results 



in the fermionic language. A full analysis of the IR 
properties of the dynamically bosonized RG flow fol- 
lows in Sec. |VII[ For the search of a fixed point, we 
hence end up with a system of six coupled nonlinear 
equations for the five couplings in the SYM regime 
(ra^, Ai, A2, h^, Nfh v /m v ) and the anomalous dimension 
t]^. The stability matrix dfii/ dgj\ g --> at a fixed point 
g* = {m*^ ,X\,\l,h*^ ,Nfh*y /m*y) has five eigenvalues, 
which we refer to as — 0i, . . . , — O5, in ascending order ac- 
cording to their real part. In the large- Nf limit the fixed- 
point structure can be mapped out analytically and we 
discover precisely the known structure from the fermionic 
flow |49) : the Thirring universality class is governed by 
a UV fixed point having one IR relevant direction with 
critical exponent Q% = 1 and Oi>2 < 0. It is located 
in the pure fermion-vector sector, i.e., Nfh*^ /mj 2 = 
and NfK*y lm*y > 0. Moreover, the fixed-point posi- 
tion exactly coincides with the Thirring fixed point in 
the fermionic RG BS1. 



h* 2 



$; = jVf-Sr = 0, 



9v 



NX 



Sir 2 



24 F) (0) : 



for Nf — > co, 



for Nf —> 00, 



(69) 



(70) 



for the dimensionless renormalized fermionic couplings 
g^ = Z~ 2 h d ~ 2 g$ and gv = Z~ 2 k d ~ 2 gv ■ For finite Nf we 
discover deviations from the fermionic fixed-point struc- 
ture. In Fig. [4] results from the bosonized description are 
shown as black lines, and the corresponding fermionic re- 
sults from [3S] are plotted as gray lines. These deviations 
are not unexpected, since the fermionic description differs 
from the partially bosonized description in two respects: 
first, the fluctuation-induced kinetic terms of the bosons 
describe momentum-dependencies of the fermionic four- 
point functions. These momentum-dependencies are ne- 
glected in a fermionic derivative expansion used in [49] . 
In addition to this advantage of the partially bosonized 
description, there is also a disadvantage arising from the 
fact that the fluctuations are no longer Fierz-complete 
in the partially bosonized description even in the point- 
like limit, as will be demonstrated explicitly in the next 
section. 

The question which description is quantitatively more 
accurate cannot be answered a priori. There are, how- 
ever, a couple of indications that the fermionic descrip- 
tion is more reliable at least in the UV: First, as al- 
ready observed in [49], the universal aspects of the UV 
fixed-point structure in the fermionic language do not 
depend on the choice of regulator, as it should be. If 
momentum-dependencies were important, we would ex- 
pect a strong scheme dependence of the universal observ- 
ables on approximate RG calculations. This UV stability 
of fermionic flows has been observed in a variety of con- 
texts [10011101] . Second, the fixed-point properties in the 
fermionic and partially bosonized description differ both 
quantitatively and qualitatively, as is visible in Fig. [4] for 
the Thirring fixed point. Also, in the bosonized language, 
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FIG. 4. Left panel: nonuniversal UV fixed-point values of the Thirring fixed point for the bosonic (black) and fermionic RG 
flow (gray, compiled from [49]) for the linear regulator. In the bosonic language already for JVf > 4 the fixed point is located 
in the pure fermion- vector sector, while in the fermionic language this is the case only for large JVf. Fermionic and bosonic 
descriptions coincide only for Nt — > oo but approach each other again for small N{. Right panel: largest critical exponent Bi 
and subleading exponent O2 (inset). 



we find two further interacting fixed points with two or 
more relevant directions (while there is only one fixed 
point with two relevant directions in the fermionic de- 
scription). In fact, we observe a qualitative breakdown 
of the partially bosonized description for intermediate 
Nf-. for N{ < 4 the additional fixed points are located in 
the pure fermion-scalar sector h v /m v = and the pure 
fermion- vector sector h^/rn^ = 0, respectively. Since 
d t h^ oc 0(hV) and d t h v <x 0(h v ) both sectors are in 
fact invariant under RG transformations. 

The Thirring fixed point, which for JVf < 4 has compo- 
nents in both sectors, h^/m^ > and h v /m v > 0, hits 
the fixed point in the pure fermion-vector sector once 
Nf — > 4. The fixed-point position as a function of ATf 
is therefore nonanalytic at Nf = 4. For Nf > 4 the 
Thirring fixed-point stays in the fermion-vector subspace 
h^/rrA = 0, coinciding with the fermionic-RG fixed point 
only for Nf 00. In Fig. [4] we have plotted the critical 
exponents for the RG relevant (Oi) and the RG irrele- 
vant direction (62), showing good agreement in the for- 
mer case while in the latter case, in particular for small 
Nf, large deviations occur. 

To summarize: for large JVf > 1, we find a qualita- 
tive behavior which is consistent between the fermionic 
and the partially bosonized description. The scalar sec- 
tor decouples for Nf ^> 1 and only the fermion-vector 



interactions Nfhy/r 



-g matter. For smaller JVf, 



the quantitative deviations between both descriptions 
increase and particularly the fixed-point structure be- 
comes rather different, culminating at a nonanalyticity 
at JVf = 4. We interpret this behavior as a breakdown 
of the partially bosonized description arising from the 
Fierz incompleteness as shown below. Incidentally, the 
agreement between both descriptions improves again for 
smaller JVf, where the flow in both cases is dominated 
by a strong fermion-scalar coupling Nfh^/m 2 ^ ~ — 2g. 
This single-channel dominance appears to alleviate the 
problem of Fierz incompleteness. This observation also 



lends support to the recent single-channel analysis of the 
JVf = 1 system [33]. 

We conclude that the "Fierz ambiguity" is a se- 
vere problem for partially bosonized formulations of the 
Thirring model precisely in the region where different 
channels compete with each other. For the location of the 
quantum phase transition in terms of a critical fermion 
number, a solution of this problem is mandatory. For ap- 
plications of the Thirring model at a fixed given fermion 
number, say JVf = 2 for graphene, a solution of this prob- 
lem may still exert a sizable quantitative influence. A so- 
lution in terms of a dynamically bosonized RG flow will 
be presented in the next section. 



VII. DYNAMICALLY BOSONIZED RG FLOW 



Although the bosonic partition function (25) is fully 
equivalent to the original fermion theory (21), the cor- 



responding leading-order truncations of the effective ac- 
tion are not. The reason is that the four-fermi couplings, 
though absent in the bare action due to the Hubbard- 
Stratonovich trick, are again generated by the box di- 
agrams displayed in Fig. [5] However, an inclusion of 
g c j >l k and gv,k in the truncation (37) does not seem very 



appealing, as it introduces a redundancy in the effec- 
tive action. Instead, we employ dynamical bosonization 
[67l [7311741 1102] as it has already been successfully used to 
resolve mean-field ambiguities |66] or implement particle- 
hole fluctuations in ultracold fermi gases [103] , Here we 
use the simplified approaclj^] proposed in [67] . The idea 
is to perform a Hubbard-Stratonovich transformation at 
each RG step, such that all newly generated four-fermi 



2 This approach follows from an approximation to an exact equa- 
tion [73| where the neglected terms are parametrically suppressed 
for the present application I74| . 
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FIG. 5. Box diagrams contributing to the flow of g$ and gv ■ Solid lines are fermions, dashed lines are scalar fields, and wiggly 
lines vector fields. Doubled inner lines denote full propagators G k = {Pko + Rk)' 1 ■ 



interactions are again reexpressed in terms of the bosonic 
interactions. In this way, the four-fermi couplings g^ t k 
and gv,k vanish at all scales. The bosonic fields then 
necessarily become scale dependent. We use the follow- 



ing field redefinitions 

4> l Ldk = C - Kx'x 1 )^^, (f> A = <j>, (71) 

V^.t-dk = V^ k + (tfcr^Suv^ V A = V, (72) 

with to be determined functions uwyfc. Note that we 
keep the fermion fields fixed. For scale-dependent bosonic 
fields the flow equation for the effective average action is 
modified, 



d k T k [^V k ] = d k T k [^V k ]\, k y k + / E? |pW? + / ^I^^W 

= \ STr + i / ^%^(x^)d k ^, k - f ^^{t^)^ (73) 



where the first term is evaluated for fixed fields and hence 
leads to the standard flow of with <j> a and V\ replaced 
by (j) k and V k , respectively [67]. We have suppressed the 
additional dependence of T k on the fermion fields \ and 
X for brevity. Projecting onto the boson couplings, we 
find the beta functions 



d t u = d t u\ 



4>k;Vk 



gV 2 = d t v 2 



d t m 2 v = d t mv\^ Vu , d t hl = dth%\^ Vi + u'dtu^k, 



<t>k,v h ' 

<kk,Vk 



(74) 



d t h v = d t h 



v 1 4>k, Vk 



(75) 
m v d t uj Vikl 
(76) 



i.e., the scale-dependent bosonization changes only the 
flow of hi and h v and leaves the other beta functions in 
the bosonic sector invariant. For the four-fermi couplings 



we obtain 



dt94, = d t94>\ 
d t gv = d t gv 



4>k,V k 



4>k,V k 



h^dtiOfjy^, 
- hydtCOv.k 



Choosing 



d t UJ V ,, 



(77) 
(78) 



(79) 



where ^ := &t9<l>\^ htVh , P gv := d t gv\ tj)k y k , establishes 
that gv. k and g^ k vanish at all scales, if absent at the 
UV scale k = A. The beta functions f3 94> and (3 gv are 
straightforwardly obtained by suitable projections of the 
box diagrams in Fig. [5j 
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+ i i ™ )d (%h%u' + 2Ku";r 1x 



4v d a 



2Nf 2 £ { ™ )d (^h 2 ,u'; Vx ,^) + (27V f 2 - l) ^(fih^u' + ^u f ; Vx , m ) 

h\ + iv d afj v ( [ i® )d {f- { hli m\ + vk\ r) x ,ri v )hy 
2N { 2 e™ d (^hl,u',m v]Vx ,^) + (2N f 2 - l) l [ ^ )d {^-h%u' + f-u f ,m 2 v ;fj x ,^) 

+ 4 F iT )d (j%^ u ' + 2 ™"> m h h l h v 



(3) 



<t>/V 



(80) 



with to be determined (possibly iVf-dependent) constants 

(i) 

a^jy In the %SB regime the fermions couple also 
to the expectation value k and there are more terms 
~ nl^fi f ( K tij,/Nt, . . . ) h^yh 2 ^. They are suppressed for 
small k <§; 1 and bounded from above for large n ^> 1 
since Kl ( ™ )d {. . . ) ~ «/(l + nh^/Nt) ■ e^ )d {. ■■)■ We 
expect that such terms do not take a significant influ- 
ence on the flow. For simplicity, we will omit them in the 

following. Instead of evaluating the diagrams in Fig. IB] 

(i) . 

explicitly, we can determine the u^/v ky taking advan- 
tage of the fact [551 EZ] that the dynamically bosonized 
flow in the pointlike limit Zv t k ~~ ^ an d Z^^ — ¥ exactly 
coincides with the fermionic flow computed in [49], 









d t 


Nt 










d t 













= -2dtg\ 



= d t (g-g)\ 



<Pk,Vk' 



(81) 
(82) 



This fixes 



0/V 



uniquely and establishes an exact map- 
ping of the fermionic fixed-point structure |49j onto the 
bosonized language in the pointlike limit. Beyond the 
pointlike approximation, the bosonized RG permits to 
reliably run toward and into the xSB regime, allowing us 
to predict the desired IR values of, for instance, fermion 
mass or order parameter as a function of Nf. 



VIII. UV STRUCTURE AND FIXED POINTS 

A. Large-iVf limit 

In the limit of infinite flavor number Nt — > oo, the 
flow equations simplify considerably. For the dynami- 
cally bosonized flow in the symmetric regime we have 



(83) 



dtiriy 



(-2 + Vv )m 2 v + —£^(0; Vx ) (Nfh 2 v ) 



d t (NtXi) = (-l + 2^)(JV f Ai)- 



— f 



(P) 



(-l + 2^)X 2 --^i ( p(0;r lx )hl, 



d t {N it i) = (-l + 2 Vv )(N lf i) 



d t (Ntis) = (• 



dth 2 = (- 



15tt 2 
3tt 2 2 



4 F) (0; % ) (N f h 2 v 



(85) 
(86) 

(87) 



f T] V ) (N { u) 
(0;r, x ) (N { h 2 y) 




(0,m v ;j] x ,j] V ) (Nfhy) 



d t (Nth 



V4> 
Vv 



(-1 

2 
3^2 

4 

3^2 
0, 



2g x ) (N f h 2 v ) , 
m ( p(0;r, x ) (Nth 2 ,) , 



(89) 
(90) 

(91) 

(92) 
(93) 



where we have multiplied the flow equations with suitable 
factors of Nf in order to simplify the large- Nf counting 
of orders. Hence, in this limit the fermion wave func- 
tion renormalization does not flow, r\ x — —dt In Z x ^ = 0, 
whereas the flow of the bosonic wave function renormal- 
izations is nonzero. 

Let us search for fixed points: assuming a fixed-point 
value m*? ^ 0, we are forced to conclude from Eq. (83) 



that the scalar anomalous dimension at the fixed point 
obeys rf^ — 2. Similarly, Eq. (90) requires for an inter- 



(84) 



acting fixed point with hy ^ for the vector anomalous 
dimension rfy = 1. In fact, this is the large- Nf result 
found in [20] • With these values of the anomalous dimen- 
sions, the full set of fixed-point values can be computed 
analytically, see Table [TTJ 

At this point, it may be worthwhile to make a few 
comments. For taking the limit Nf — > oo, we have im- 
plicitly assumed a specific iVf-scaling of the couplings, 
defining the corresponding 't Hooft couplings. Since all 
rescaled couplings acquire finite (in particular nonzero) 
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TABLE II. Left columns: nonuniversal fixed-point couplings for the linear regulator and various flavor numbers Nf. Right 
columns: universal correlation length exponent v, subleading exponent ui, and anomalous dimension Rough error estimates 
arise from the comparison to sharp cutoff results. For iVf > 6 we do not expect chiral symmetry breaking due to vector 
channel domination, such that the critical exponents of the Thirring fixed point do not have the same physical meaning for 
the long-range physics. This is indicated by the gray font of the corresponding critical exponents; see Sec. |IX| For Nf = 2, a 
meaningful estimate of the systematic error for v is not available due to sharp-cutoff artifacts, see text. 





m%INf* 




A 2 


h? 


N { h&/m& 


V 






2 


0.132 


5.85 


24.68 


15.96 


1.99 


2.4(?) 


1.0(2) 


1.4(7) 


3 


0.265 


16.94 


39.84 


19.36 


6.11 


1.22(2) 


1.5(2) 


1.6(3) 


4 


0.271 


20.95 


44.45 


20.61 


8.83 


1.084(5) 


1.7(1) 


1.7(2) 


5 


0.266 


22.70 


46.55 


21.41 


10.83 


1.043(8) 


1.74(7) 


1.8(1) 


6 


0.259 


23.56 


47.70 


21.95 


12.35 


1.026(6) 


1.80(4) 


1.85(8) 


8 


0.247 


24.34 


48.86 


22.59 


14.44 


1.012(3) 


1.84(2) 


1.91(5) 


10 


0.238 


24.68 


49.42 


22.93 


15.79 


1.007(2) 


1.79(4) 


1.94(3) 


12 


0.231 


24.85 


49.73 


23.14 


16.74 


1.004(9) 


1.72(4) 


1.95(2) 


25 


0.209 


25.16 


50.33 


23.54 


19.39 


1.0006(1) 


1.423(1) 


1.988(2) 


100 


0.187 


25.26 


50.52 


23.68 


21.45 


1.0000(1) 


1.126(1) 


1.999(1) 


oo 


8 

45 

~ 0.178 


64-7^ 
25 

~ 25.27 


128^ 
25 

~ 50.53 


~ 23^69 


9tv' a 

~ 22.21 


112 



fixed-point values, the corresponding rescalings are self- 
consistent and the large- Nf limit is well-defined. For the 
original four-fermi couplings = Nfh^/m 2 ^ and gy = 
Nfhy /rriy this implies the large- Nf scaling g$ = 0(1/ Nf) 
and gy = 0(1), already known from the fermionic 
flow |35]. In fact, the fixed-point positions of and 
gv for the dynamically bosonized flow in the large-A^f 
limit and fermionic flow exactly coincide, provided the 
same regulator is employed. We note, however, that 
the equivalence in general does not hold for the flow be- 
yond the fixed-point regime, since kinetic terms of the 
bosonic fields are generated during the flow according to 
their finite anomalous dimensions. For instance, the large 
anomalous dimension rj*^ = 2 indicates a rapid flow of the 
scalar kinetic term. 

Consider the second contribution to the beta function 
of the scalar- fermion coupling /i| in Eq. ( 89 ) , being pro- 
portional to m|. This is exactly the contribution from 
the diagrams in Fig. [5] i.e., the large- Nf reminiscence of 
the scale-dependent Hubbard-Stratonovich transforma- 
tion. Without this contribution, the large- Nf fixed-point 
equations cannot be solved with finite couplings. In other 
words, dynamical bosonization is crucial to find the cor- 
rect large- Nf behavior in the scalar sector in the present 
model. We have already seen an indication for this in 
Sec.|VIt where the scalar sector completely decoupled for 
Nf > 4 using only standard partial bosonization. 

For the given set of integer anomalous dimensions, 
the large- Nf fixed-point values for the vector-vector cou- 
pling (i and the vector-scalar coupling v are in fact nega- 
tive. We attribute this to the fact that we neglected the 
momentum-dependent contribution cx C,k in Eq. (37 1 as 
well as higher order terms in the effective potential. In 
any case, in the large- Nf limit fi and v do not feed back 
into the flow of the remaining couplings. For consistency, 
we again evaluate in what follows the flow equations for 



a pointlikc current-current interaction, Zy,k ~^ 0, and 
neglect the vector-vector selfinteractions fi and vector- 
scalar interactions v, analogous to Sec. |VI| This is ex- 
pected to be a reasonable approximation as long as the 
vector mass my does not become too small. 



B. Fixed-point structure at general Nf 

Beyond the large- Nf limit, we evaluate the fixed-point 
equations numerically, both for linear and sharp cutoff. 
Again, we recover the known UV structure: there is one 
interacting Thirring fixed point for all 1 < Nf < oo, 
having only one IR relevant direction. [Recall that the 
case Nf — 1 has been explicitly excluded in the deriva- 
tion of our bosonic flow equations, cf. Eq. ([44]).] For 
small Nf the Thirring fixed point is located close to the 
pure scalar channel subspacc, whereas for Nf — > oo the 
scalar-fermion coupling becomes negligible. As we have 
also observed in the partially bosonized formulation in 
Sec. |VI[ the fixed-point position in this limit again ex- 
actly coincides with the fermionic fixed point. For small 
Nf, we find deviations from the fermionic UV structure. 
As the dynamically bosonized flow and the fermionic one 
have an identical fixed-point structure in the point-like 
limit (by construction), these deviations are fully related 
to the momentum dependence in the bosonized scalar 
channel. The latter is particularly important once the 
scalar coupling becomes long range, i.e., for small m?. 
In other words, the bosons dynamically become fluctuat- 
ing relevant degrees of freedom. We depict the fixed- 
point positions for the linear regulator in Fig. [6] to- 
gether with the corresponding values for the fermionic 
flow compiled from [49 for comparison. Explicit values 
are given in Table [TT] If the fixed point corresponds to 
a second-order phase transition as we expect it for small 
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FIG. 6. Left panel: nonuniversal UV fixed-point values of dynamically bosonized RG flow for the linear regulator. The couplings 
have been multiplied with suitable Nf factors; for better visibility, the scalar mass has been multiplied with an additional factor 
50. Right panel: Comparison of dynamically bosonized RG (black) with fermionic RG (gray). The improvement due to 
dynamical bosonization becomes obvious by comparing this plot with the left panel of Fig. [4] 
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FIG. 7. Left panel: correlation length exponent v extracted from the linearized flow in the fixed-point regime with linear (solid) 
and sharp regulator (dashed). The difference (gray shaded band) serves as rough error estimate. The critical exponent is to a 
large extent independent of the regulator, while the uncertainty increases for Nf \ 2, as shown in the inset. For increasing Nf 
it approaches rapidly the fermionic result l/0i = 1. Right panel: corrections-to-scaling exponent w. 



Nf, the critical behavior is uniquely determined by the 
fixed-point regime, where the flow can be linearized. Let 
g = {rr&, Ai, A 2 , h^, Nfh v /m v ) denote the vector of our 
(generalized) couplings. The scaling of the correlation 
length in the vicinity of the critical point is 

£ = m^ R « \S 9 r(l + b ± \6 9 r + ■■■). (94) 

with 8g :— g\ — g C r measuring the distance from critical- 
ity ( "reduced temperature" ) and the renormalized scalar 
mass R = lim^o m^k 2 . By again denoting the small- 
est eigenvalue (being negative) of the stability matrix 
d(3i/dgj\ g * with —Q\ < and the second smallest eigen- 
value (being positive) with —02 > we have for the 
correlation length exponent and the first subleading ex- 
ponent |104j 

z/ = l/6i>0, and uj = -6 2 > 0. (95) 

At the critical point Sg — 0, where the correlation length 
diverges, the asymptotic behavior of the scalar two-point 



function is determined by the anomalous dimension r\*^ = 
V4>{9*) as 

(<p(xM0)) oc ^ d \ +vy (96) 

The critical exponents extracted from the flow in the 
fixed-point regime are shown in Figs. [7] and [8] We ex- 
pect the values, as listed in Table [TTJ obtained with the 
linear regulator to represent our most accurate results 
and use the difference to the sharp-cutoff results as a 
rough estimate on the truncation-induced error. We find 
that the correlation-length exponent is to a large extent 
regulator- independent, Av/v < 1 . . . 2% for Nf > 2.1. 
Near Nf = 2, the sharp-cutoff results deviate from the 
linear regulator by up to a factor of two. We inter- 
pret this behavior as a large artifact of the sharp cut- 
off and consider the result using the linear regulator as 
our best estimate. For this particular case of Nf = 2, 
it is however difficult to estimate the systematic error. 
As expected, slightly larger but controlled deviations 
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FIG. 8. Anomalous dimension t]^ at the fixed point, with 
linear regulator (solid) and sharp regulator (dashed). 




- 0.6 - 0.4 - 0.2 

In(fc/fc*) 

FIG. 9. RG evolution of renormalized masses in xSB regime 
for k < k* and Nf = 2. We stop the flow when the radial 
mass w,r jP approaches its maximum (vertical line). 



between linear and sharp cutoff occur for the anoma- 
lous dimension and the corrections-to-scaling exponent, 
A^/fj^, Aw/w < 10... 15%. Explicit values are again 
given in Table [TTJ. This hierarchy of accuracy is well- 
known from RG studies of scalar models based on the 
derivative expansion [701 USB 1105] . 



IX. IR BEHAVIOR AND iV f -CONTROLLED 
QUANTUM PHASE TRANSITION 

A. Critical fermion number 

If we start the flow for small Nf with initial UV cou- 
plings close to the fixed point, we find that the scalar 
mass eventually vanishes at some scale k* , indicating the 
spontaneous breakdown of chiral symmetry. In the fol- 
lowing, we will refer to k* as "xSB scale" . Continuing the 
flow for k < k* in the xSB regime the fermions become 
massive with renormalized mass rh 2 ^ f = Nf^ 1 nh 2 k 2 and 
the scalar sector consists of one radial mode with renor- 
malized mass m\ = 2/iAifc 2 , 2Nf 2 — 1 massive modes 
with m\ r = Nf nXzh 2 , and 2Nf 2 massless Goldstone 
modes. In the deep IR k — > 0, we expect all massive 
modes as well as eventually the Goldstone modes to de- 
couple, leading to IR predictions for the mass spectrum. 

However, our truncation is not able to resolve this de- 
coupling, as all couplings among the different modes are 
approximated by the same expansion coefficient of the 
effective potential. In fact, we naively find that the di- 
mensionless parameters run into an attractive IR fixed 
point, such that the dimensionful masses run into a max- 
imum and eventually decrease again for small k -C k*, 
see Fig. [9| This effect is a well-known artifact of our 
treatment of the effective potential, as it artificially cou- 
ples the Goldstone modes to the flow of all other scalar 
operators. This problem can technically be solved by an 
adapted choice of field coordinates jl06j . Here, we fol- 
low a more immediate strategy, and simply stop the flow 
before it enters the artificial IR fixed-point regime. For 



our quantitative estimates, we stop the flow at the max- 
imum of the radial mass fh\ . We have checked that 
the following results to a large extent do not depend on 
this choice. We note that the critical behavior in terms 
of the exponents v and rf^ remains, of course, unaffected 
by this, since it is solely determined by the UV structure 
of the theory. 

Once the physical scale has been set, for instance, by 
measuring the value of the radial mass m| , we can com- 
pare the dynamically generated masses among different 
Nf. We have already shown the renormalized fermion 
mass m\ f in units of the radial mass m\ p in Fig. |2j In 



Fig. 10 we now depict the renormalized scalar mass mj T 



in the same units. For increasing flavor number we ob- 
serve that both fh 2 ^ f and decrease (apart from some 
interesting nonmonotonic behavior of the fermion mass 
for Nf ~ 3) and eventually vanish for some critical flavor 
number Nf /*• N"' . This is a clear and expected signa- 
ture for a quantum phase transition governed by flavor 
number Nf. 

we find that the vector-fermion 



Directly at N f = N? 

,2 /™2 



f 

coupling Nfhy/rriy diverges at the same time as the 



scalar mass mi k approaches zero, i.e., at the x^B scale 
k* . We expect that this divergence is an artifact of our 
truncation which will be stabilized by higher terms in 
the vector sector [for instance, the terms oc flk, v^, f& 
in Eq. (37)]. Here, we interpret the divergence as an 
indication for "vector-boson-dominance" in the IR, in- 
hibiting xSB. If the divergence was real, the model could 
exhibit dynamical Lorentz symmetry breaking [91H95] . 
For larger Nf > iV f cr we observe the divergence already 
at higher scales k > k* , inhibiting the flow to enter the 
xSB regime. 

The value of iV f cr can thus be obtained in two ways: 
First, we observe for increasing Nf a sharp decrease of 
the logarithm of the fermion mass once Nf /~ N". Sec- 
ond, we look for the largest Nf below which the vector- 
fermion coupling Nfhy/niy remains finite for all scales 
in the interval k g [fc*,A]. Above this Nf, the vector- 
fermion coupling diverges before the scalar mass m ( p_k 
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FIG. 10. Scalar mass ra R r in units of the radial mass fh K p on 
logarithmic scale (linear cutoff). Gray shaded area: estimates 
for critical flavor number iV f cr ~ 5.1 (linear cutoff) and 5.8 
(sharp cutoff). 



FIG. 11. RG evolution of dimensionless couplings for the 
"pure" Thirring model with Nf = 2 and initial couplings 
K,a = 0, h^ A = 0, and two different gv,A = 

9c 



10° 



"4,a 

Nfhy A /my A just above (black) and below (gray) g CI . Near 
ln(fc/A) ~ —5 the flow approaches the UV fixed-point regime. 
The black solid line corresponds to above and k below the 
XSB scale ln(fe*/A) — —15, respectively. 



approaches zero. For the linear cutoff we find N" ~ 5.07 
with the former method and JV f cr ~ 5.1 with the lat- 
ter, agreeing with each other on the level of numerical 
precision. For the sharp cutoff, our truncation unfortu- 
nately does not allow to enter the xSB regime: we find 
Xi.k* < at the xSB scale, indicating the requirement 
of higher terms in the polynomial expansion of the ef- 
fective potential ( 43 1 . This is in line with optimization 
considerations [731 186] : for a nonoptimized regulator (as 
is the sharp cutoff) a higher expansion order is needed 
to achieve similar predictive power. For the scalar and 
fermion sectors, the linear regulator by contrast satisfies 
standard optimization criteria [75|l86j. Nonetheless, also 
in the sharp-cutoff scheme, we can determine the flavor 
number above which the vector-fermion coupling diverges 
before m^k ~^ 0- From this criterion, we find N" ~ 5.8. 
Identifying the cutoff dependence with a rough error es- 
timate yields our result for the critical flavor number 



N™ ~ 5.1(7). 



(97) 



This prediction represents one of the central results of 
this work. 



B. Thirring universality class 

In the fermionic RG 49], we have shown that the 
"pure" Thirring model, defined by a microscopic action 
including only current-current interaction, is in the uni- 
versality class defined by the interacting Thirring fixed 
point with one IR relevant coupling. Of course, this state- 
ment ultimately only holds for our specific regularization 
scheme, since fixed-point positions itself are nonuniversal. 
In the bosonized formulation developed here we can fur- 
ther explore the universality class of the "pure" Thirring 



model by starting the flow with initial couplings 

Z 1 ° cm 0A> 1 ' Ai, A <l = 1,2), ^ A <1, 

(98) 

i.e., with a decoupled scalar sector. In our explicit com- 
putations we use Za, a = 10~ 6 , rru, a = 10 6 , = 0, 



h 2 ^ A = 0. However, universality guarantees that our re- 
sults are independent of the exact values for the initial 
couplings as long as the system is governed by the corre- 
sponding fixed point. With the functional RG flow, we 
can explicitly verify universality. 



Fig. 



11 shows the RG 
2 for 



evolution of the dimensionless couplings for Nf 
two different initial values of the vector-fermion coupling 
gv = Nfhy/rriy just above and below the critical value 
g cr . Though we start with a decoupled scalar sector, the 
scalar couplings are rapidly generated by the RG flow, as 
already known from the purely fermionic RG. For initial 
couplings close to criticality, the flow runs from the pure 
Thirring axis into the fixed-point regime, which is exactly 
given by the Thirring fixed-point couplings in Table [XXJ 
cf. Fig. [TT] Thus, the critical behavior of a convention- 
ally defined "pure" Thirring model without scalar chan- 
nel in the bare action is indeed given by the universality 
class of our Thirring fixed point. For completeness, we 
show in Fig. [12] the critical coupling g cr for which the 
flow approaches the Thirring fixed point as a function of 
Nf, to be compared with previous results in Ref. [49] {^] 
As in the fermionic calculation, we do not find a sharp 
decrease of l/g cr as Nf /*• N". This observation is in 
agreement with lattice simulations [2TH3T] but disagrees 



3 We use the opportunity to point out that the values for l/g C r 
plotted in |49j include an erroneous factor of 2. This has already 
been corrected in 1501. 
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FIG. 12. Critical coupling l/g CI for the "pure" Thirring model 
with initially decoupled scalar sector as a function of N{. For 
Nf > 5.1 (gray area) we find no signature of chiral symmetry 
breaking (dashed line). 



fixed-point regime (see Table |TT| , together with the hy- 
perscaling assumption: 

|(d-2 + ry;)=2.771, i/(2 - t£) = 1.539, ^ = 2.360, 

(102) 



in excellent agreement with Eqs. (101). However, hyper- 



scaling is to a large extent also fulfilled for the convention- 
ally defined "pure" Thirring model, cf. Eqs. ( |100[ ) with 
(102 1. We attribute the small (but significant) hyper- 



scaling violations to the presence of other fixed points: 
in particular the critical behavior could be influenced 
by fixed point B, which (at least in the pointlike limit, 
see Fig. [IJ is located in the broader vicinity of the pure 
Thirring axis and potentially describes a different univer- 
sality class. We will further elaborate on this aspect in 
the following section, see also Fig. [Ml 



with earlier analytical estimates [24] [25] . A rapid varia- 
tion of 1/ <? cr would have been expected if the many-flavor 
quantum phase transition was induced by a change in the 
UV fixed-point structure. This, however, is not the case 
in the Thirring model where the transition is rather due 
to a competition between the vector and scalar channel. 
Since the Thirring fixed point is present for all Nf, we can 
still determine a "would-be" critical coupling g cr even for 



N t > Nf 1 ' (dashed line in Fig. [12 1 



In order to establish the connection between the nega- 
tive eigenvalues of the stability matrix 0^ and the critical 
exponents in Eq. ( 95 ) we have assumed that the propa- 



gator has a scaling form, from which one infers the hy- 
perscaling relations [104] 



(99) 



The hyperscaling assumption can in fact explicitly be 
checked by computing directly the critical behavior of 
the order parameter (ip) oc |^g|' 3 , inverse susceptibility 
(unrenormalized mass) x _1 = "ip <x l^ffl 7 ) an( i inverse 
correlation length (renormalized mass) £ = oc 
\Sg\ u . As a function of the distance to criticality, we 
indeed find that the expected linear behavior on a double- 
log plot is excellently fulfilled; see Fig. [13] for the case 
of the "pure" Thirring model (initially decoupled scalar 
sector) and Nf = 2. The slopes of the regression lines 
for initial couplings on the pure Thirring axis are (for 
iV f = 2) 



P = 2.769, 7 = 1.553, 



2.364, 



(100) 



whereas if we start the flow directly near the Thirring 
fixed point, we obtain 



P = 2.771, 7 = 1.539, v = 2.361. 



(101) 



The values should be compared with the critical expo- 
nents obtained from the linearized flow in the Thirring 



X. COMPARISON WITH PREVIOUS STUDIES 

The diversity of results for the 3d Thirring model ob- 
tained so far in the literature calls for a careful compar- 
ison of our findings with those derived with other tech- 
niques. With hindsight, this variety of partly contradic- 
tory results arises from the somewhat unexpected com- 
plexity of the model. We have identified several sources 
for this complexity: (i) an involved UV fixed point struc- 
ture (three instead of just one UV fixed point with the 
Thirring fixed point being off the pure Thirring axis) , (ii) 
a complex bound state spectrum giving rise to effective 
vector and scalar degrees of freedom, (iii) a symmetry 
breaking mechanism with competing vector and scalar 
channels. All issues need to be taken care of in order to 
arrive at conclusive answers. 

Let us start with a comparison with the large- Nf ex- 
pansion. With standard large- Nf techniques, renormaliz- 
ability of the "pure" 3d Thirring model with only current- 
current interaction oc gvi^Pl^) 2 (i- e -> without scalar 
channel) has been shown to hold if and only if a reg- 
ularization scheme is employed in which the vector field 
propagator remains purely transverse on the quantum 
level, defining an interacting UV fixed point 1!) 22"] 1107] . 
In Ref . [5D] it is found to leading order in 1 /Nf (in terms 
of our notation) 



dtgv = 9v 1 



9v 
9*v 



W = ^, Vx=0(l/N { ), 
9v 

(103) 



which is exactly the large-iVf behavior of our flow 
equations (83 )— (93 (: if the fixed point corresponds 



to a second-order phase transition, the corresponding 
correlation- length exponent is v = 1/0 = 1 where 







-dp/dg 



and the vector field anomalous dimen- 



sion is i] V (9v) = L 

The seeming resemblance to gauge theories has been 
taken even more seriously: Yang [107] claims that the 
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FIG. 13. Critical behavior of "pure" Thirring model with N{ = 2: order parameter (ip) oc \5gf (left panel) and inverse 
susceptibility x _1 — TO? oc |<5<?| 7 (right panel) in double-log plot, showing the expected power law. The slope of the regression 
line (gray dashed) is /3 = 2.769 and 7 = 1.553, respectively. Left inset: for very large coupling on the pure Thirring axis the 
order parameter decreases again, see Sec. [X] 



partially bosonized Thirring model is equivalent to (a 
gauge-fixed version of) a U(l) gauge theory, in which 
the mass of the vector boson is generated by the Higgs 
mechanism. Insisting on the gauge symmetry, it is there- 
with argued that the coupling gv cannot be renormal- 
ized. Consequently, the beta function would be vanishing 
for any value of gv [2H EE] ■ 

From our functional RG viewpoint, nonpcrturbative 
renormalizability does neither rely on a specific regu- 
larization scheme nor a resemblance or equivalence to 
a gauge theory. All that is needed is a non-Gaufiian UV 
fixed point (or a line of fixed points) with suitable proper- 
ties to render the system asymptotically safe [5TJ[55]. We 
observe such a fixed point both in the purely fermionic 
description which does not have a gauge symmetry as 
well as in the bosonized language where the vector mass 
term breaks gauge symmetry manifestly. Still, it is inter- 
esting to see that our large- Nf flow equations ( 83 )— (93 1, 
resemble that of a gauge-fixed theory: the vector mass 
does not feed back into the vector sector and the second 
kinetic term (~ Ay j.) looks like a gauge-fixing term that 
is even locked to the kinetic term in our truncation. 

Whether or not the RG flow is eventually attracted 
by a gauge-invariant theory in the long-range limit still 
remains an open question. In our approach, this could 
be investigated by measuring the flow with respect to 
the hyperplane of actions in theory space that sat- 
isfy the (regulator-modified) Ward-Takahashi identity 
[731 [73J 1108] . We believe, however, that this question 
is unrelated to that of chiral symmetry breaking, as the 
latter occurs in the chiral scalar sector which does not 
participate in the local symmetry. 

The formulation of the Thirring model as a gauge the- 
ory is also at the basis of the analytical studies [2"2U2"5"] . 
which attempt to solve the Dyson-Schwinger equations 
(DSE) by setting the full vector propagator to its large- 
Nf form, and the full vertex to the bare vertex (known 
as ladder approximation). All such studies observe the 



existence of a nontrivial solution corresponding to chiral 
symmetry breaking for small values of Nf < N". How- 
ever, their predictions for the value of iV f cr and the critical 
behavior of fcrmion mass and chiral order parameter dif- 
fer significantly: Ref. E2] reports -/V f cr ~ 3.24. Hong and 
Park [23 claim that symmetry breaking should persist 
for any value of Nf, i.e., Nf T = 00. In Ref. [24] an essen- 
tially singular behavior close to Nf 1 ~ 4.32 is found, 



mjt^f oc A exp 



2tt 



for Nf < N 



r ■ 



(104) 



that is to say, a phase transition of infinite order. Such 
an essential scaling law is known from a number of gauge 
theories, such as quenched QED 4 [42 , 109] 1110) and 
QED 3 with Nf fermion flavors [BHH1 QH EES], as well as 
many-flavor QCD 4 [4"2l |4U1 [TntiTTB] . In this context, the 
scaling law is often referred to as Miransky scaling; the 
scenario has also been termed (pseudo-)conformal phase 
transition [32] ■ Essential scaling has also been found 
in two-dimensional statistical systems, such as the XY 
model, where the scenario is called Kosterlitz-Thouless 
phase transition [114H117] . As has been pointed out re- 
cently [79] 1118j , the general mechanism being responsible 
for essential scaling is the merger of two RG fixed points 
|100| I119[ I120j and their subsequent disappearance into 
the complex plane, or a running of a fixed point off to 
zero or infinite coupling. By contrast, we find that the 
UV Thirring fixed point persists for any Nf, hence there 
is no basis for and consequently no observation of essen- 
tial scaling in our work. 

By constructing an effective potential for the chiral 
order parameter (ip) oc (V^A), to leading order in 1/Nf, 
Kondo |26j reports a second-order phase transition, 



A f cr 
~Nf~ 



for Nf < 7V f cr , 



(105) 



21 



where for d = 3 he finds the critical exponent 6=1 and 
iV f cr = 2. We have been motivated by these considera- 
tions to attempt a fit to our results for order parameter 
and fermion mass. In Fig. [3j we depict (tp) /mR.,p versus 
the combination (N" /Nf) — 1 in a double- log plot. In 
fact, the behavior is very well compatible with a power- 
law scaling corresponding to a second-order phase tran- 
sition at N™. However, the linear fit (gray line) gives 
b ~ 0.44 in qualitative agreement with the Kondo sce- 
nario p6| but a quantitatively differing estimate for the 
exponent. We emphasize that for a quantum phase tran- 
sition at N" occurring due to competing condensation 
channels, the corresponding critical exponents do not 
have to coincide with those determining the phase tran- 
sition with the coupling as control parameter. In partic- 
ular, there is no reason that /3\ Ncr and b should coincide. 
In any case, a fit to an essential scaling behavior analo- 
gous to (104) both for order parameter and fermion mass 



is much less successful and is not supported by our re- 
sults. 

Chiral symmetry breaking in the "pure" 3d Thirring 
model has also been investigated on the lattice [27U32] . 
One extensive study [3T] finds N™ = 6.6(1), the bare 
value of which we consider as still consistent with our 
result, see below. Moreover, the study also analyzes the 
critical behavior close to the phase transition, both for 
fixed Nf < N" as a function of the bare coupling g and 
for fixed g > g cr as a function of Nf. In the latter case, 
the data is fitted to an equation of state of the form 



= A 



(Nf- N^+CL-n (W)P+B(iPi)) 



\p+j 



(106) 



with m being the bare fermion mass (explicit symmetry- 
breaking parameter) , L the linear extent of the system, 
and A, B, C some constants. Close to criticality in the 
chiral (m — > 0) and thermodynamic (L — > oo) limit, 
Eq. (106) in fact reduces to the Kondo-formula (105). 
It is interesting to notice that the fit reported in [31] 
yields b ~ 0.37, which could be considered as roughly 
compatible with our result b ~ 0.44. 

However, striking discrepancies occur, when one com- 
pares the critical behavior for fixed Nf < N" with the 
bare coupling as control parameter. In particular, in 
the sequence of studies [2"8"H3"T] it is found that the ex- 
ponent S lies between 8 « 2.8 (Nf — 2) and 8 rs 5.8 
(Nf — 6), i.e., most notably, S increases with Nf. The 
result for Nf = 2 has been confirmed in a recent work [35] 
based on a new algorithmic approach, where the same 
lattice action has been employed directly in the massless- 
fermion limit with manifest U(1)®U(1) chiral symmetry 
The resulting anomalous dimension reads rf^ = 0.65(1), 
which (under the hyperscaling assumption) is equivalent 
to 8 = (5 - r/p/(l + rfy) = 2.64(2). By contrast, we find 
(cf. Table [H]) 8 « 1.5 (Nf = 2) and 6 « 1.1 (Nf = 6), 
that is to say, S decreases with Nf and the values lie well 
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FIG. 14. Critical exponents 8 and v for a theory defined at 
the fixed point B with more than one RG relevant direction; 
to be compared with Fig. 5 of Ref. |31j . 

below the above-mentioned lattice predictions^] 

The lattice simulations rely on a microscopic defini- 
tion of the Thirring model which is fixed by a bare ac- 
tion including only the Thirring-like current-current in- 
teraction. It is therefore a priori not clear, whether the 
critical behavior in the lattice models is in fact given by 
the Thirring fixed point, cf. the discussion in Ref. |49) . 
Our analysis in the preceding section |VIII| supports the 
implicit assumption that bare actions on the Thirring 
axis belong to the Thirring universality class for all rele- 
vant Nf (cf. Fig. 11). Nonetheless, since fixed-point po- 
sitions are scheme-dependent, our results may not be di- 
rectly transferable to the lattice theory. In particular, 
it is known from the flow in the pointlike limit, that 
a microscopic formulation, being fixed on the Thirring 
axis, could, for instance, be influenced by the fixed point 
B, which has more than one RG relevant direction (see 
Fig. [I]). In order to get a glimpse of how such a situation 
could change the corresponding critical behavior we com- 
pute in the dynamically bosonized formulation the criti- 
cal exponents associated with fixed point B. We empha- 
size however that a theory defined at B in any case cannot 
be fixed by just one bare parameter. Naively using the 
hyperscaling relation^ (6) = [5 - rfy(B)]/[l + rf^B)] we 
can relate the critical exponent 8(B) to the anomalous 
dimension rf^(B) at B, see Fig. 
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The values for 8(B) 
lie well above those for the Thirring fixed point. Most 
notably, 8(B) in fact now increases with iVf, in qualita- 
tive consistency with the behavior shown in Fig. 5 of |31j . 
For comparison, we also depict the correlation length ex- 
ponent v = l/0i, associated with the largest negative 



4 As a side remark we note that the DSE studies 24, 25 point to 
5 = 2 (Nf < Nf) and 8 = 1 (Nf = Nf). (This fact has been 
debated in Ref. 1311 .) The numerical similarity to our results is 
probably accidental, since the nature of the transition reported in 
these studies substantially differs from ours; see Eq. | |104} above. 
Of course, hyperscaling may not hold at a fixed point with two 
relevant directions. 



22 



eigenvalue of the stability matrix, i.e., the strongest RG 
relevant direction. In order to test this scenario lattice 
simulations with actions containing both invariant oper- 
ators ~ gv as well as are needed. 

However, we would like to propose yet another way 
to resolve this seeming contradiction between our results 
and lattice studies: All quoted lattice studies implement 
staggered fermions, in which case the chiral limit is more 
easily accessible. Exactly the same lattice action em- 
ployed in the Thirring studies has also been used in sim- 
ulations of the chiral U(l) ® U(l) Gross- Neveu (xGN) 
model in three dimensions |121) . yielding v = 0.88(8) 
and = 0.46(11), probably consistent with the Nf = 2 
lattice Thirring results v = 0.85(1), r\\ = 0.65(1) [32] 
and v = 0.71(4), r\\ = 0.60(2) [25]. In fact, the symme- 
try breaking pattern of the lattice Thirring model in the 
staggered fermion formulation is |29j 



2.75(9) for N st3 



1 



[We emphasize how- 



U(iV s tagg) ® U(7V st agg) -> U(JV s tagg), 



(107) 



where N stagg is the number of staggered fermion flavors. 
It is related to the number of continuum four-component 
fermions by Nf = 2N stagg [122J. It is an open question 
whether the breaking pattern of the continuum Thirring 
model, 



U(2iV £ ) -> U(iV £ ) ® U(iVf), 



(108) 



is restored in the continuum limit of the lattice Thirring 
model; see [2H1 [29] for a discussion. This scenario 
could be checked by a careful analysis of the low-energy 
spectrum in the broken phase: the number of Gold- 



stone modes for the breaking pattern (1071 would be 
NiL™ = N f 2 /4, whereas for (fl08|) one expects 2A f 2 



stagg ly t 

massless modes. Alternatively, one could examine the 
spectrum of the Dirac operator on the lattice: a neces- 
sary condition for the restoration of the full U(2Af ) sym- 
metry in the continuum limit is that the lattice Dirac 
operator has to exhibit the same four-fold degeneracy 
of the continuum Dirac operator. In fact, for a similar 
strongly-coupled (2+l)-dimcnsional model (in the con- 
text of graphene), such a comparison reveals significant 
breaking of the U(2iVf ) symmetry, even in the vicinity of 
the second-order phase transition point [123] . Unfortu- 
nately, we do not know of any such study for the Thirring 
model so far. If the continuum Thirring pattern ( 108 1 is 



not restored in the simulations, then the lattice Thirring 
result could possibly describe the xGN universality class 
with breaking pattern ( 107| . In other words, the results 
for critical exponents would differ, simply because the 
continuum and the lattice models would not be in the 
same universality class. 

We have been motivated by these consideration to com- 
pare the lattice results also to results for a fermionic 
U(A L ) (gi U(A R ) model [SOI W\- For N L = N R = 1, 
a functional RG study has provided the result <5 x gn = 
2.58(9). Interestingly enough, we observe that these 
xGN findings are in fact well consistent with the gen- 
eral trend reported for the lattice Thirring model, e.g., 



ever that N stagg = 1 corresponds to two four-component 
continuum flavors, whereas the U(l) <£> U(l) model in 
Refs. [SOI [M] is defined with solely one four-component 
fermion.] In Refs. [2IB 1121) the question has been raised 
whether for Nf — 2 the distinction between xGN and 
Thirring model might be unimportant and both mod- 
els might actually coincide. In such a scenario the fixed 
point pertaining to the x>3B phase transition would inci- 
dentally lie in a U(4)-symmetric IR attractive subspace 
of the theory space. While certainly possible, we know of 
no argument supporting such a conjecture; it is also not 
favored by our results for the xGN and Thirring model. 

In order to clarify these questions, it would be very in- 
teresting to simulate the Thirring model in a formulation 
in which either the restoration of the U(2iVf) symmetry 
can be explicitly verified in the continuum limit or in 
which it is manifest. In the two-component fermion for- 



mulation [Eq. (pi l] the latter option might be pursued us- 
ing standard methods (e.g., Wilson fermions |124j ). since 
there is no notion of chirality in this formulation and 
the U(2ATf) is just a flavor symmetry. Alternatively, in 
the four-component language [Eq. ( 14 )] , one could em- 
ploy more advanced chiral-fermion techniques, such as 
domain-wall [125H127] or overlap |128j fermions, both of 
which are explicitly chirally symmetric and deserve to 
be studied on their own rights due to their relevance to 
lattice QCD 4 . 

We would finally like to comment on the unortho- 
dox behavior of the chiral condensate far from the crit- 
ical point, as reported in the lattice study [3T]. There, 
it has been found that ('tpip) decreases again for large 
(current-current) coupling gy. The corresponding peak 
position at gy — gy,max is independent of both lattice 
volume and bare fermion mass, indicating that its ori- 
gin is at the UV scale. The effect has then been at- 
tributed to the fact that the vector propagator in the 
lattice regularization violates the transversity condition, 
the latter being crucial to the renormalizability of the 
1/iVf expansion. The impact of the extra divergence 
can however be absorbed by a coupling renormalization 
g' v = gv(l — gv / gv.Hm) 1 , such that the strong coupling 
limit is recovered at gy = gy,n m . gv,\im is thereupon 
identified with the peak position <?y,max' 

By contrast, our RG results suggest that the 3d 
Thirring model can in fact be renormalized nonpertur- 
batively without insisting on a transverse vector prop- 
agator, if one allows for a microscopic definition in the 
two-dimensional coupling plane spanned, for instance, by 
the couplings gy and g^. If this remains robust beyond 
our approximation, it could provide a natural explana- 
tion for the nonmonotonic behavior of the condensate as 
a function of gy: Since the numerical value of any IR ob- 
servable (in units of a fixed UV cutoff A) is first and fore- 
most given by the "duration" of RG time t = ln(fc/A) be- 
fore the flow freezes out or enters the IR regime, it could 
well be that for a theory defined on the pure Thirring 
axis with bare g$ — an IR quantity decreases again 
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for large gv far from criticality. Loosely speaking, just 
because we start the flow with a large current-current in- 
teraction does not necessarily mean that we are closer to 
the xSB regime. In order to check these considerations 
we have computed the order parameter for different bare 
couplings gv on the Thirring axis and an initially decou- 
pled scalar sector, showing indeed a maximum far from 
criticality gv ^> 3cr and a subsequent decrease; see inset 

of Fig. m 



XI. CONCLUSIONS 

We have studied chiral symmetry breaking in the 3d 
Thirring model with Nf fermion flavors. Using the func- 
tional RG equation for the effective average action, we 
have investigated the RG flow of the system parametrized 
in terms of the fundamental fermionic as well as com- 
posite bosonic degrees of freedom. We have analyzed 
both UV structure and IR behavior of the theory in the 
usual formulation with fixed fields as well as for dynami- 
cally bosonized fields, that is to say, by applying a scale- 
dependent Hubbard-Stratonovich transformation. Both 
formulations show the existence of a Thirring fixed point 
with one RG relevant direction. For small Nf it is located 
close to the scalar-channel subspace, while it approaches 
the vector-channel subspace (the pure Thirring axis) for 
increasing Nf. All this is well compatible with previous 
findings or indications in the purely fermionic descrip- 
tion [49] . 

For entering the symmetry-broken regime, however, 
the description using composite bosonic fields is highly 
advantageous. Whereas partial bosonization has be- 
come a standard tool for functional RG analyses of 
strongly-correlated fermion systems, a quantitatively re- 
liable study of the 3d Thirring model appears to require 
dynamical bosonization [671 1731 1741 1102] . The reason is 
that the physics of competing channels is technically af- 
fected by the Fierz ambiguity representing a challenge for 
many techniques based on Hubbard-Stratonovich trans- 
formations. The 3d Thirring model therefore serves as 
a paradigm example for the competing-channel problem 
and its resolution through dynamical bosonization. 

For the present model, we have confirmed previous re- 
sults from DSE studies and Monte-Carlo simulations in- 
dicating that chiral symmetry breaking, corresponding 
to a condensation in the scalar channel, occurs for large 
enough coupling if the number of fermion flavors Nf is 
smaller than a critical value iV f cr , while it is shown to 
be absent for Nf > iV f cr , independently of the coupling. 
We obtain the estimate N" ~ 5.1(7), where we have used 
regulator dependencies as an indicator for the systematic 
error of our approximation. 

Since all previous studies we know of rely on a mi- 
croscopic definition of the 3d Thirring model with pure 
current-current (Thirring-like) interaction, the mecha- 
nism behind the quantum phase transition at N" has 
so far been rather unclear. The present work has shown 



for the first time that the critical flavor number may oc- 
cur due to a competition between different condensation 
channels: the NJL-type scalar channel on the one hand, 
dominating for small Nf and triggering xSB, and the 
Thirring- type vector channel on the other hand, dom- 
inating for large Nf and inhibiting xSH. In particular, 
our results exclude a quantum phase transition mecha- 
nism based on fixed-point annihilation that would im- 
ply essential scaling of order parameters as a function of 



iV f c 



N { . 



The RG approach is particularly useful for predicting 
the critical phenomena associated with a given univer- 
sality class. We have computed the critical behavior for 
fixed Nf < N" as a function of the bare coupling in 
terms of the exponents u, 77^, /3, and 7 both for a micro- 
scopic definition of the model at the Thirring fixed point 
as well as on the pure Thirring axis ("pure" Thirring 
model). For Nf = 2, we have verified that the hyper- 
scaling relations are fulfilled for UV complete models de- 
fined at the Thirring fixed point. Models with initial 
conditions on the pure Thirring axis necessarily inherit a 
certain degree of nonuniversality which becomes visible 
in small violations of hyperscaling, presumably induced 
by the vicinity of another fixed point (fixed point B in our 
notation). Strictly speaking, the model starting on the 
pure Thirring axis is not itself a UV complete quantum 
field theory. We can, however, think of it as belonging 
to an RG trajectory (line of constant physics) that em- 
anates from fixed point B and thus has a definite UV 
completion. For the theory defined at the Thirring fixed 
point, we have also computed the corrections-to-scaling 
exponent oj. 

Close to the quantum critical point Nf /* iV f cr , we 
have discussed the scaling behavior as a function of Nf 
(for fixed couplings above their critical values). For the 
quantum phase transition as a function of N" — Nf, we 
have computed the exponent b ~ 0.44 ("magnetization 
exponent"), implying that there is clear evidence for the 
transition to be of second order. We have also computed 
the dynamically generated fermion mass f (fermion 
gap) and the r-mode mass fh^ in units of the radial 
mass Trip (inverse healing length) as a function of Nf < 
iV f cr . In principle, these are inherent predictions of our 
analysis and could be verified by correlator measurements 
in lattice simulations. Since our IR analysis is affected 
by an artificial non-decoupling of Goldstone modes, our 
quantitative estimate should, however, not be taken too 
literally. We consider our results on critical exponents as 
our most accurate quantitative predictions. 

At the present stage, the long-range behavior in the 
vector-dominated phase for large Nf > Nf 1 is difficult 
to resolve. Our approximation with pointlike vector 
channel does not allow a reliable prediction of the IR 
properties of the system for Nf > iV f cr , since for small 
vector mass m v < 0(1) momentum-dependent terms 
cx Zv,k)Av,k)Gk m the effective action (37) become im- 



portant as we can read off from the vector anomalous 
dimension growing large. In this sense, our estimate for 
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the critical flavor number might be a lower bound to the 
true value of iV f cr . A larger critical flavor number could, 
for instance, arise from a nontrivial interplay of dynami- 
cal vector and scalar channels driving the system to crit- 
icality also for larger flavor number. Thus, we consider 
our findings for TV™ to be still compatible with the lattice 
results [31] . pointing to jV f cr ~ 6.6. In any case however, 
it would be very interesting to see in how far the inclu- 
sion of momentum-dependent terms in the vector sector 
would modify the IR flow of the large- N{ theory. This is 
in particular true for the term oc Cfe^u^tAK/) which has 
no analogue in scalar field theories and thus could lead 
to a qualitatively different IR behavior as compared to 
the latter. Furthermore, at this point we also cannot ex- 
clude with certainty that this term does not change the 
UV structure of our theory in the large- Nf limit. This 
deserves further investigation. 

Whereas most [22" 1 124H3T ] (but not all [23]) of the pre- 
vious studies agree at least on the very existence of a 
critical flavor number at order N" ~ 0(2 ... 7), the na- 
ture of the phase transition has so far been a substantially 
delicate issue. Based on our detailed predictions, in par- 
ticular for the transition at fixed N{ < Nf 1 as a function 
of the coupling, we believe that it may now be possible 
to resolve the discrepancies in the literature. For that 
purpose, we propose an independent investigation of the 
critical behavior of the 3d Thirring model, for instance, 
by means of a Monte Carlo simulation with a lattice ac- 
tion exhibiting manifest U(2iVf) symmetry. 
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Appendix: Loop integrals 

For self-containedness, we give here the explicit results 
for the loop integrals for linear and sharp cutoff. The 



regulator functions Rk present in the Wetterich equa- 
tion ( 23 ) may be written in terms of dimensionless shape 



functions r k via 



R<t>,k{q) = Z^^r^^q 2 ), (A.l) 
Rip,k{q) = -Z^,. k dr^,. k (q 2 ), (A. 2) 



with collective bosonic and fermionic fields <fi = 
(cj) ab , Vp, . . . ) and -0 = (ip a , ■ ■ • ), respectively. The lin- 
ear cutoff, which satisfies an optimization criterion 129 , 
is defined as 



':( q ')= 0(k 2 -q 2 ) 




i e^-g*). 



(A.3) 
(A.4) 



The sharp cutoff is defined as the a — > oo limit of the 
class of regulators given by 



(A.5) 



t'i — )-i)e(k 2 - q 2 ), 




(A.6) 



where we demand for definiteness that the sharp-cutoff 
limit a — y oo is to be taken after the integration over the 
internal momentum q, i.e., after the substitution into the 
threshold functions 11301. 



The one-loop structure of the Wetterich equation guar- 
antees that the flow equations can always be written in 
terms of single integrals — the threshold functions, which 
encode the details of the regularization scheme. Their 
definitions are 



£^ /F)d (uj;^ N ,) = -k- d d t J dsn:*- 1 log [P^x) + cok 2 ] , 

) = ^^rf^V^/vO = -\k 2n - d ~d t £ dxxi- 1 [P H ^{x)+uk 2 



ffiffai.wa;^,^) = ~-k 2 ^ +n ^- d d t J^ dxxi- 1 [P^x) + Wl * 



21 -™i 



Pv(x) + w 2 /s 



,21 -"2 



4™£(wi.W2;»?0,'7*) = --k 2 ^+ n ^- d d t J o dxxl- 1 [P^+^k 2 ]^ 1 [P^x) +uj 2 k 2 ] 



21 -«2 



(A.7) 
(A.8) 
(A.9) 
(A.10) 
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■ 21 ~™2 



1 f 

x [P4x) + u 3 k 2 f n3 , (A.11) 



TO 2 B 2 d ( w 1^2;770) 



A: 6 - d 5t / da; a; 5 



x 2 



1 /"°° 
mf )d {uj-^) = --k & - d d t dx 

2 Jo 

1 Z" 00 
m { p d (Lo;r,^) = --k 4 - d d t dx 
1 Jo 



X 2 



+ 1 



'P (x)+ Wl fc 2 

1 

: P^(x)+wfc 2 _ 
1 + r^(x) 



P^x) + u 2 k 2 



P^(jc)+wfc 2 



,4 l + r^(s) 



P/,(a;) +wifc 2 x P 0/y (x) + w 2 fc 2 ' 



(A.12) 
(A.13) 
(A.14) 
(A.15) 



with m € N and where we have suppressed the scale index fc for the sake of simplicity. Moreover, we have abbreviated 
the (inverse) regularized propagator parts by 



P<t,{x) := x [1 + 7>(x)] , 
For the integrations, we have substituted q 2 i->- x, viz., 



P/,(x) := x [1 + 



JjSf= 4Vd j = 2Vd j _1 



with v d := iVol(S" i " 1 )/(27r)' i = [2 d + 1 n d / 2 T(d/2)]- 1 . d t is defined to act only on the regulator's t-dependence, 



dt ■= V / <b- 



5P^(x>) 



(A.16) 
(A.17) 

(A.18) 
(A.19) 



Both the linear and sharp regulators have the very convenient feature that all loop integrals can be performed explicitly. 
For the linear cutoff the results are 



£ o B)d (w;770) 



4 B) V;»fo) 



i 



i - 



d + 2 1+uj 



ri4> 



4 BB „) 2 d K, W2 ;^,r, v ) 

m^l d {u 1 ,u 2 -,r]<t,) 
mf )d {u;ri^) 



d + 2) (l + cj)"+ r 
?7V 



rf + iy (i+w) 



n+l ' 



ni 



d + 2 / 1 + lu 1 



+ 



n 2 



d+l) 



d + 2 J 1 + LU 2 

n 2 



d + 2) 1 + LJ 2 



(1 +Wi)"i(l +w 2 )™ 2 ' 
1 

(1 +W 1 )"l(l +LU 2 ) n2 ' 



1-il 



d + 1 J 1 + wi 



+ 1- 



d + 2y V: 



™2 



+ 



"3 \ 



d + 2yv i + w 2 i + w 3 j 

X 



f 



(i +wi)"i(i +w 2 )" 2 (i +u 3 y 



(1+ Wl ) 2 (l+ W2 ) 2 ' 

1 



(1+0,)' 



(A.20) 
(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 

(A.26) 
(A.27) 
(A.28) 
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1 



m ( p d (uj; rfy) 

(FB)d, \ A »fa/V \ 

mi 2 {cji,u 2 ;V^,V^/v) = 1 1- I 
For the sharp cutoff we find 



1 - W 



(1+w) 4 d-2 (1+cj) 3 V 2d -4 4/ 



1 



\2' 



(1+ Wl )(l+ W2 )2- 



if /¥)d {io- = - log(l + w) + ^ B/F)d (0; 



e /F)d (-;^) 

4^(wi,W2;^,f^) 

(B)d, \ 
m 2 F)d (cj;^) 

(F)d, \ 



1 + w)" ' 



1 



1 +Wi)"i(l + cj 2 )™ 2 ' 
1 

1 + + cj 2 )™ 2 ' 

1 



1 +CJl)"l(l +CJ 2 )" 2 (1 +CJ 3 )™3 



1 + W1 ) 2 (1 + W 2 ) 2 ' 
1 



1 



1+w) 4 ' 

1 



1 + W1 )(1 + W 2 ) 2 ' 



(A.29) 
(A.30) 

(A.31) 
(A.32) 

(A.33) 

(A.34) 

(A.35) 

(A.36) 

(A.37) 

(A.38) 

(A.39) 
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